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ABSTRACT 

This  paper  presents  the  results  of  the  development  of 
graphic  simulations  of  stochastic  models  using  a  micro- 
computer. Included  are  four  computer  programs  simulating 
situations  which  are  based  on  an  underlying  Poisson  process. 
The  programs  are  written  in  BASIC  programming  language  for  an 
Apple    II    Plus   computer   with    a   minimum   of    48K   of   memory. 
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I.  INTRODUCTION 

This  paper  presents  the  results  of  an  effort  to  develop 
visual  simulations  of  a  variety  of  stochastic  models  via  a 
relatively  new  media,  the  microcomputer. 

The  intent  of  the  effort  was  to  create  a  real  time 
perception  of  the  models  in  action;  to  give  the  viewer  an 
experience  that  could  not  otherwise  be  attained  through  text 
and  classroom  study.  The  microcomputer  is  suitable  for  this 
project  due  to  its  graphic  capability,  portability  and 
accessabili ty . 

The  theme  for  this  project  is  the  Poisson  process.  Four 
computer  programs  are  provided  that  simulate  the  process 
either  by  itself  or  in  applications.  Represented  are  a 
simple  and  a  superimposed  Poisson  process,  an  M/M/1  queue,  a 
dam  model  and  a  damage  model. 

Descriptions  of  the  programs  and  instructions  to  operate 
them  are  found  in  the  User  Guide.  The  remainder  of  this 
paper  will  touch  upon  several  aspects  of  the  project  and 
programs  that  may  be  of  interest  to  the  general  reader  or  of 
importance  to  those  who  may  wish  to  pursue  work  of  this 
nature.   Program  listings  are  found  in  Appendix  A, 


II.   DISCUSSION 

A.  THE  COMPUTER 

An  Apple  II  Plus  computer  with  48K  of  memory  was  used  for 
this  project.  An  Apple  was  selected  because  of  its 
relatively  easy  accessabil  i  ty  to  anyone  who  may  wish  to  run 
the  programs.  In  addition,  it  is  expected  that  the  programs 
of  this  package  will  run  with  little  or  no  modification  on 
future  generations  of  Apple  hardware.  The  programs  will  run 
on  an  Apple  III  in  the  emulation  mode. 

B.  THE  LANGUAGE 

The  programs  in  this  package  are  written  in  Applesoft 
BASIC,  a  subset  of  Microsoft  BASIC  with  floating  point  and 
special  graphics  features.  The  selection  of  this  language 
over  other  languages  available,  (including  PASCAL,  and 
FORTRAN)  ,      was      a      matter      of      convenience.  Apple     is     a 

'BASIC-speaking'  computer,  and  there  is  a  great  deal  of 
commercial  software  support  for  developing  programs,  graphics 
and    other    screen   displays. 

Although  BASIC  was  used,  Pascal  was  the  first  choice  for 
this      project.  Its      structured      nature      provides     high 

legibility,  aiding  the  interested  reader  as  well  as  the 
original  programmer  in  understanding  the  program.  The  author 
switched    to    BASIC    after    discovering    that   many   specialized 


graphics    features    required    for    the    programs   were    not   avail- 
able    in    commercial    Pascal    graphics    packages.      Rather    than 
make      a      compromise     on      the     graphic     presentation     of     the 
simulations,    the    author   opted    for    an   alternate    language. 

The  type  of  programming  language  used  is  not  critical. 
Speed  is  not  essential,  nor  are  there  any  special 
capabilities  needed  other  than  those  provided  by  Applesoft 
BASIC.  Most      of      the      programming      involves     use     of     the 

computers*     resident    graphics   package    and    a    second    special 
graphics    package,    to    be   described   next. 

C.  SOFTWARE  AIDS 

Two  commercial  software  packages  were  used  to  aid  in 
producing  the  programs  of  this  project;  Printer  II  by 
Computer  Systems  Design,  required  to  format  numerical  output 
in  two  of  the  programs,  and  the  'Hi-Res  Character  Generator' 
(HRCG)  found  in  the  Applesoft  Tool  Kit  produced  by  Apple 
Computer,  Inc.  The  Printer  II  package  adds  a  final  touch  to 
the  programs  and  is  not  absolutely  necessary  to  successful 
production.  The  HRCG  was  essential.  Without  it,  many  of  the 
graphics  features  would  have  been,  at  the  least,  very 
difficult    to    produce. 

D.  PROGRAMMING    METHODS 

Program  listings  are  found  in  Appendix  A.  As  will  be 
noted,      most     of      the      programming      structure     deals     with 


preparation  and  execution  of  graphics.  All  of  the  programs 
are   structured    in   a   similiar    fashion. 

A  typical  program  may  be  divided  into  several 
subroutines.  The  initialization  portion  of  the  program 
starts  at  line  number  10000.  Other  important  parts  include; 
a  main  timing  or  clock  routine  starting  at  line  1000, 
plotting  routines  starting  at  2000,  a  keyboard  command 
routine  starting  at  3000,  and  in  two  cases  a  parameter  change 
routine    starting    at    line    4000. 

The  initialization  portion  of  the  program  defines 
character  strings,  sets  constants,  and  determines  time  to  the 
first  event  of  the  program.  The  timing  routine  is  similiar 
to  a  next-event  structured  simulation.  The  program  loops 
within  this  routine,  keeping  count  of  each  pass;  the  count 
representing  time.  In  each  pass  a  check  is  also  made  to  see 
whether  it  is  time  for  an  event  to  occur.  If  so  the  program 
jumps    to    the    plotting    routine. 

The  plotting  routine  handles  all  the  particulars  needed 
to  output  on  to  the  screen  the  effect  of  an  event  occurring, 
for  example,  a  customer  arriving  at  a  queue,  or  an  automobile 
passing  by  a  particular  point  along  a  road.  After  the  event 
occurs  and  while  still  in  the  plotting  routine,  a  random 
number  is  obtained  which  is  subsequently  converted  to 
represent  time  to  the  next  event.  Then,  the  program  jumps 
back    to    the    timing    routine. 


The  keyboard  routine  is  simply  a  'control  center'  to 
filter  commands  intercepted  from  the  keyboard.  Keyboard 
commands  not  recognized  are  ignored.  The  parameter  change 
routine  starting  on  line  4000,  is  provided  to  allow  the  user 
to    change   mean    interarrival    times    in    two   of    the    programs. 

E.   GRAPHIC  SIMULATIONS 

There  are  53,760  locations  arranged  in  a  280  by  192  array 
on  the  monitor  screen  available  for  program  output.  This 
presents  a  problem  in  truly  representing  a  continuous  random 
variable  such  as  exponential  interarrival  times  of  the 
Poisson  process.  Fortunately,  for  the  purposes  of  this 
project,  this  limitation  is  minor.  The  interarrival  times  as 
represented  in  the  programs  appear  to  belong  to  a  Poisson 
process   and    the    appearance    is   considered    sufficient. 

All  of  the  programs  in  one  way  or  another  describe  events 
and  the  interarrival  times  between  events.  The  interarrival 
times  shown  or  perceived  are  in  fact  geometric;  the  discrete 
analog  of  the  exponential.  The  random  number  representing 
time  to  the  next  event  is  rounded  up  to  the  next  highest 
integer  to  do  this..  The  average  interarrival  time  is 
slightly  higher  than  the  average  that  would  result  if  the 
random  numbers  were  truly  exponential.  The  visual  effect  on 
the    screen    from    this    alteration    is   minimal. 
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F.   RAND0r4  NUMBERS  AND  TRANSFORMATIONS 

The  following  paragraphs  briefly  explain  how  the  programs 
use  random  numbers  to  create  Poisson  simulations.  The  source 
of  random  numbers  is  the  BASIC  function  RND.  RND  supplies 
uniform  (0,1)  random  numbers  which  are  then  transformed  in 
various  ways  to  suit  program  objectives.  The  sequence  of 
random  numbers  used  in  the  programs  are  not  repeatable, 
though  this  option  was  available. 

Each  of  the  four  programs  present  simulations  based  on  an 
underlying  Poisson  process,  and  therefore,  all  use 
exponential  random  numbers.  To  obtain  the  exponentials, 
uniform  (0,1)  random  numbers  are  transformed  through  the 
inverse  probability  method.  The  variable  'y'  is  set  equal  to 
the  uniform  random  number,  and  the  variable  'x*  is  assigned 
the  resultant  exponential  random  number.  The  constant  'L'  is 
the  rate  parameter; 

F(x)  =  y  =  1  -  sxp(-  L  *  X) 

1  -  y  =  exp(-  L  *  X) 

ln(l  -  y)  =  -   L  *  X 

X  =  -  ln(l  -  y)  /  L 

or;  x  =  -ln(y)/L 

The  above  transformation  is,  of  course,  valid;  F(x)  is  the 
distribution  function  of  a  continuous  random  variable  and  the 
inverse  distribution  function  exists,  therefore,  the 
probability  density  of  F{x)  is  uniform  on  the  interval  (0,1). 
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The  program  labeled  'Poisson  Process'  presents 
simulations  of  both  a  simple  and  a  superimposed  Poisson 
process.  For  the  simple  process,  exponential  random  numbers 
are  obtained  one  at  a  time.  Initially,  one  exponential 
random  number  is  generated.  As  the  program  progresses,  a 
time  counter  is  incremented  and  eventually  becomes  equivalent 
to  the  number,  whereupon,  an  event  is  simulated.  The  program 
then  generates  another  exponential  random  number,  the  counter 
is  reset  and  the  procedure  begins  again. 

In  the  superimposed  mode,  two  Poisson  processes  are 
active,  each  with  a  different  rate  parameter,  (call  them  Type 
I  and  Type  II  with  respective  parameters  LI  and  L2)  .  Two 
random  numbers  are  obtained  initially,  one  for  each  type.  As 
the  program  progresses,  a  time  counter  is  incremented  until 
the  value  of  the  counter  becomes  equal  to  the  smaller  of  the 
two  numbers.  Then,  a  simulation  of  the  type  of  event  that 
the  smaller  number  is  associated  with,  occurs.  After  the 
simulation,  this  smaller  number  is  replaced  with  a  new  random 
number.  The  next  event,  (Type  I  or  Type  II),  will  occur  when 
the  counter  becomes  equivalent  to  the  smaller  of  the  two 
numbers  now  on  hand;  the  new  one  or  the  one  that  was  not 
involved  with  the  last  simulation. 

Should  the  values  of  the  two  random  numbers  ever  be 
equal,  a  uniform  (0,1)  random  number  is  obtained  and  used  to 
determine  which  of  the  two  types  of  events  will  occur  first. 
If  the  number  is  less  than  a  predetermined  value  'p',  or\e 
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type  of  event  occur,  and  if  not,  the  second  type  of  event 
will  occur.  The  value  'p'  is  generated  by  the  following 
equation  where  'LI'  and  'L2'  are  the  rates  of  the  two  Poisson 
processes; 

p   =    LI    /    (LI    +    L2) 

The  value  'p'  is  the  probability  that  a  Type  I  event  will 
occur  before  a  Type  II  event.  Let  'T'  be  the  time  to  the 
next  event  and  'Tl'  and  'T2'  be  independent  random  variables 
corresponding  to  the  respective  times  that  the  Type  I  and 
Type  II  events  will  occur;  then,  the  time  to  the  next  event 
will  be; 


T   = 


Tl  with   probability        p 


T2  with   probability        1   -    p 


There  is  another  transformation  that  takes  place  when 
mean  interarrival  times  are  changed  by  the  user  through 
keyboard      commands.  In      that      transformation,     the     mean 

interarrival  times  displayed  on  the  screen  are  non  those  used 
in  the  transformation  equations;  they  are  converted  first 
through  an  equation  relating  the  exponential  distribution  to 
the  geometric  distribution.  Below,  the  variable  'GM'  is 
assigned  the  value  obtained  from  the  keyboard  input,  (and 
displayed    on     the     screen),     and     the    variable     'EM'     is    the 
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parameter  used  in  the  transformation  equation,  (EM  =  mean 
interarrival  time  =  1  /  L;   L  =  rate  of  the  process); 

EM  =  -  ln(l  -  1  /  GM) 

Thus  if  'GM'  is  made  equal  to  5,  then  'EM'  will  be  assigned 
the  value  4,48142012.  This  procedure  is  included  in  the 
program  only  to  synchronize  the  results  presented  in  the 
statistics  display  with  the  time  plot. 

The  M/M/1  queue  model  is  built  very  much  like  the  program 
•Poisson  Process'  in  that  two  exponential  random  numbers  are 
needed,  one  for  interarrival  times,  and  one  for  service 
times.  The  numbers  are  obtained  in  the  same  way  and  the 
exponential  to  geometric  conversion  is  identical.  There  is 
no  concern  for  the  possibility  of  a  simultaneous  occurrence 
of  an  arrival  event  and  a  departure  event.  Thus,  if  the  two 
random  variables  should  become  equal  the  time  plot  included 
in  the  program  indicates  a  simultaneous  arrival  and 
departure.  Character  animations  representing  arrivals  and 
departures  allow  for  the  arrival  to  always  occur  first. 

The  two  remaining  programs  involve  the  compound  Poisson 
process.  In  the  dam  model,  the  two  random  numbers  are  used 
whenever  an  event  occurs  and  therefore,  are  both  replaced  at 
the  same  time.  One  represents  the  interarrival  time  to  the 
next  event  (arrival  of  a  cloud)  and  forms  the  Poisson 
process.  The  second  is  an  independent  and  identically 
distributed  random  variable  independent  of  the  interarrival 
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times,  (the  amount  of  rain  that  falls  from  each  cloud) .  This 
second  random  variable  is  also  exponential.  Random  variables 
with  other  distributions  could  have  been  used  but  the 
exponential    was    suitable    for    the   graphic    presentations. 

The  second  compound  Poisson  process  is  found  in  the 
damage  model  program.  Here,  as  with  the  dam  model,  two 
random  numbers  are  obtained  at  the  same  time;  both  are 
exponential.  One  is  used  to  form  the  Poisson  process  and 
represents  times  to  successive  events  that  cause  damage.  The 
amount  of  damage  incurred  is  represented  by  the  second  random 
variable.  Again,  a  random  variable  with  some  different 
distribution   could    have   been    used    for    amount   of   damage. 

In  both  compound  Poisson  process  models  the  random 
variables  are  not  converted  from  geometric  to  exponential  or 
vice   versa. 

G.   THE  RANDOM  NUMBER  GENERATOR 

As  noted  previously,  the  BASIC  function  RND  was  used  to 
obtain  uniform  (0,1)  random  numbers  for  program  operation. 
It  is  not  known  how  closely  a  sequence  of  these  numbers 
conforms  to  a  uniform  d i s t r i bu i to n .  The  nature  of  the 
programs  does  not  require  exactness  and  therefore  precision 
was  not  of  great  concern.  Nevertheless,  out  of  interest,  two 
standard  tests  for  distribution  conformity  and  independence 
were  run  to  obtain  some  measure  of  the  generators' 
suf f ic  iency . 
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The  two  tests  run  were  the  frequency  test  and  the  serial 
test.  Each  was  run  100  times  and  for  each  of  the  100 
iterations,  1500  random  numbers  were  used.  The  results  were 
favorable.  Comparison  of  the  observed  sequences  with  the 
expected  characteristics  of  the  uniform  (0,1)  distribution, 
showed  a  3%  rejection  for  the  frequency  test  and  a  5% 
rejection      for      the      serial      test.  Using      a      chi-square 

good ness-o f - f i t  test  for  each,  the  hypothesis  that  the 
distribution  from  which  random  numbers  are  obtained  for  the 
simulations  is  uniform  (0,1),  cannot  be  rejected  at  the  95% 
level    of   confidence. 

For  work  requiring  precision,  further  tests  of  RND  would, 
of  course,  have  to  be  conducted,  and  no  definite  conclusions 
concerning    its    adequacy   can   be    drawn   here. 
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A.  INTRODUCTION 

The  package  of  programs  collectively  labeled  'POISSON 
PROCESSES,'  is  designed  for  use  on  the  Apple  II  or  Apple  II 
Plus  microcomputer.  The  programs  themselves  are  designed  to 
provide  the  viewer  with  a  real  time  perception  of  Poisson 
processes  in  action.   Each  simulates  a  different  situation; 

-  'Poisson  Process,'  is  a  simple  demonstration  of  the 
occurrence  of  events  in  a  Poisson  manner, 

-  'MMIQ,'  is  the  simulation  of  a  simple  single-server 
queue  with  exponential  interarrival  times  and  exponential 
service  times  (the  M/M/1  queue). 

-  'Dam  Model,'  simulates  one  of  the  classic  applications 
of  the  compound  Poisson  process. 

-  'Damage  Model,'  represents  the  compound  Poisson  process 
applied  to  a  cumulative  damage  model. 

The  following  paragraphs  provide  an  explanation  of 
equipment  requirements,  procedures  to  follow  to  properly 
operate  the  microcomputer,  and  a  more  detailed  description  of 
each  program. 

B.  EQUIPMENT  REQUIREMENTS 

An  Apple  II  Plus  computer  (or  Apple  II  with  Applesoft 
card  or  Language  System) ,  with  at  least  48k  of  memory,  a  disk 
drive  and  controller,  and  a  monitor  (or  television)  are 
required  to  run  these  programs.   The  minimum  memory  size  is 


21 


needed  to  provide  two  'pages'  of  memory  for  graphics.  An 
Apple  is  required  not  only  because  the  programs  are  written 
in  Applesoft  Basic  but  because  the  programs  use  a  special 
high  resolution  character  generator  program  made  for  the 
Apple  by  Apple  Computer  Inc.,  to  produce  special  effects  on 
the  screen. 

C.   STARTING  UP 

To  run  the  programs,  place  the  disk  into  the  disk  drive, 
and  turn  on  the  computer  and  monitor.   The  introductory 
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Figure  1:  First  Display  of  Indroductory  Program 


program,  (FISHPOND  LOGO),  will  be  automatically  loaded  and 
run.  It  will  take  a  few  seconds  before  anything  appears  on 
the  screen,  as  the  startup  program  must  first  load  the  high 


resolution  character  generator  and  character  sets.   The  first 
presentation  will  look  like  Figure  1,  above. 

By  pressing  the  space  bar  as  directed  in  the  first 
display  the  presentation  will  change  to  that  shown  in  Figure 
2; 
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Figure  2:  Second  Display  of  Introductory  Program 

The  'Acknowledgements'  section  of  the  introductory 
program  sights  the  use  of  commercial  software  used  to  run  the 
programs.  The  'Introduction'  is  a  condensed  version  of  the 
introduction  in  this  guide  and,  therefore,  will  not  be 
discussed  here.  Through  the  'Menu,'  the  user  can  access  each 
of  the  four  programs  which  make  up  the  package. 

Once  into  a  program,  any  of  the  other  programs  may  be 
accessed  either  through  the  Menu  or  directly.  To  obtain  the 
Menu  at  any  time,  simply  press  the  'M'  key.   To  go  to  another 
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specific  program  directly,  use  the  same  key  as  called  for  in 
the  Menu. 

In  addition,  each  program  has  its  own  'Help'  program 
which  can  be  accessed  by  pressing  the  'H'  key.  These  Help 
programs  give  condensed  versions  of  the  following  program 
descriptions  as  well  as  brief  listings  of  keyboard  commands 
which  can  be  used  in  the  respective  programs. 

D.   PROGRAM  DESCRIPTIONS 

Each  program  is  described  separately  on  the  subsequent 
pages.  Each  description  includes  a  brief  explanation  of  what 
is  happening  in  the  program  and  includes  a  summary  of  the 
commands  that  can  be  used  during  the  program  operation. 

Program  1 
POISSON  PROCESS 
This  program  illustrates  a  simple  Poisson  process  in 
action.  The  graphical  representation  is  composed  of  three 
parts;  a  time  plot,  animated  characters  which  move  across  the 
screen  at  random  intervals  of  time,  and  an  optional 
statistics  display.  The  advance  of  time  is  represented  by  a 
plot  of  'tic'  (or  time)  marks.  To  better  associate  real  time 
with  program  time  consider  the  unit  of  time  as  a  'second.' 
Different  types  of  marks  are  plotted  for  one  'second',  five 
'seconds,'  and  sixty  'seconds.'  The  occurrence  of  an  event 
is  marked  with  a  diamond  shape  and  an  animated  character  (bus 
or  automobile),  moves  across  the  screen  from  left  to  right. 
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When  the  program  is  first  run  the  viewer  will  be  looking 
at  a  Poisson  process  which  is  in  fact  a  composition  of  two 
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Figure  3:  'Poisson  Process'  Screen  Display 

Poisson  processes;  one  process  for  the  transit  of  automobiles 
and  a  second  process  for  the  transit  of  busses.  This 
superimposed  process  may  be  decomposed  by  pressing  the  'T' 
key.  Subsequently,  two  lines  of  time  marks  will  appear,  the 
upper  line  for  the  auto  and  the  lower  line  for  the  bus, 
(re-compose  by  pressing  the  'U'  key). 

There  are  several  other  features  to  this  program  which 
the  viewer  may  find  useful.  First,  the  user  may  alter  the 
program  so  that  only  one  type  of  event  will  occur;  automobile 
or  bus,  (press  the  'A'  key  for  auto  only;  the  'B'  key  for  bus 
only).  Secondly,  the  user  may  change  the  mean  interarrival 
times  of  either  the  auto,  bus  or  both,  (press  the  'C  key  and 
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follow  the  instructions).  A  statistics  display  is  also 
provided,  (press  the  'D'  key  to  turn  on).  The  display  shows 
a  time  count,  and  the  theoretical  and  empirical  mean 
interarrival  times  for  the  automobile  and  bus.  The 
theoretical  means  are  built  into  the  program  and  may  be 
changed  as  noted  above.  The  empirical  means  change 
automatically  with  the  passage  of  time,  (the  respective  means 
are  identified  in  the  display  by  subscripts;  'A'  stands  for 
automobile  and  'B'  stands  for  bus;  Press  the  '0'  key  to  turn 
the  display  off  )  . 

There  are  several  ways  to  speed  up  the  program  so  that, 
for  example,  one  may  watch  the  values  of  the  empirical  means 
converge  to  those  of  the  theoretical  means.  First,  the 
program  has  two  built  in  speeds;  slow  and  fast.  The  slow 
speed  is  the  default,  (press  'F'  for  the  fast  speed  and  'S' 
to  get  back  to  the  slow  speed) .  Secondly,  the  program  will 
run  a  little  faster  without  the  vehicles  running  across  the 
screen,  (press  'N')/  and  faster  still  if  the  time  marks  do 
not  display,  (press  'N'  again;  The  statistics  display  must 
be  on  to  turn  the  time  marks  off;  Press  the  'Y'  key  twice  to 
return  all  screen  representations) . 

As  an  example,  a  key-press  sequence  might  proceed  as 
follows;  Press  'D'  for  the  statistics  display;  Press  'N'  to 
stop  the  vehicles  from  running  across  the  screen;  Press  'N' 
again  to  stop  the  'tic'  marks  from  being  plotted.  Now,  only 
the  statistics  display  will  be  showing,  (Again,  note  that  the 
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statistics  display  must  be  showing  before  the  'tic'  marks  can 
be  deleted  from  the  screen).  To  bring  everything  back, 
reverse  the  sequence;  Press  'Y',  then  'Y'  again,  then  '0'  to 
delete  the  display. 

Program  2 
MMIQ 
This  program  provides  a  visual  display  of  the  M/iM/1 
queue;  a  s  i  ng  1  e- se  rver  queue  with  exponential  interarrival 
times  and  exponential  service  times.   As  with  the  'Poisson 
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Figure  4:   'MMIQ'  Screen  Display 

Process'  model,  this  program  too  provides  a  three  part 
graphical  presentation;  animated  figures  that  move  across  the 
screen,  'tic'  marks  that  record  the  passage  of  time  and  an 
optional  statistics  display. 
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The  keyboard  commands  are  nearly  the  same.  The  user  may 
alter  the  mean  interarrival  time  and  service  time,  (press  key 
'C')/  change  the  program  speed,  ('S'  for  slow  and  'F'  for 
fast),  view  a  statistics  display,  ('D'  to  turn  on  and  '0'  to 
turn  off)  ,  and  stop  the  character  animation  and  time  plot  to 
speed  up  the  program. 

The  screen  presentation  will  show  two  lines  of  'tic' 
marks.  The  upper  line  provides  a  representation  of  time  and 
customer  arrivals,  which  are  marked  by  an  open  diamond  shape. 
The  lower  line  depicts  service  times  and  a  departure  event 
marked  by  a  solid  diamond  shape.  When  there  are  no  customers 
at  the  'service  counter'  no  service  time  marks  will  be 
plotted.  It  will  be  noted  that  the  customer  being  serviced 
is  close  to  the  'service  counter'  and  not  part  of  the  queue. 

The  statistics  display  shows  a  time  count,  the  queue  size 
(Q#)  ,  and  the  average  queue  size  (Avg  Q#)  .  The  theoretical 
means  are  not  shown.  They  may  be  determined  however,  by 
entering  the  parameter  change  routine,  (press  key  'C'). 

Program  3 

Dam  Model 
This  program  portrays  one  of  the  classic  applications  of 
the  compound  Poisson  process.  In  the  model,  clouds  arrive  in 
a  random  fashion  and  with  each  appearance  a  random  quantity 
of  rain  falls.  The  cloud  arrivals  form  a  Poisson  process, 
and  the  amounts  of  rain  from  the  clouds  are  independent  and 
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identically  distributed  random  variables,  independent  of  the 
cloud  arrival  times,  (the  rainfall  amounts  are  exponential  in 
this  model) .  There  is  an  outflow  of  water  from  the  dam  at  a 
constant  rate. 


•'•','  'i' •''','  'i' 
1 1  ' '  'I  'I  I  « '  'I  ' 


i 


Figure  5:  'Dam  Model'  Screen  Display 

There  are  few  keyboard  commands  for  this  program.  As 
with  the  two  previous  models,  there  is  a  fast  clock  time  and 
a  slow  clock  time,  ('F'  and  'S').  In  addition,  the  viewer 
may  interrupt  the  program  to  raise  or  lower  the  water  behind 
the  dam,  ('U'  to  raise  the  level,  'D'  to  lower  the  level,  and 
'N'  to  return  to  normal  program  operation).  No  statistics 
display  is  available,  nor  is  there  a  parameter  change 
routine. 
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Program  4 
DAMAGE  MODEL 

This  program  depicts  the  compound  Poisson  process  applied 
as  a  cumulative  damage  model.  The  graphical  representation 
is  composed  of  two  parts;  a  time  plot  represented  by 
horizontal  line  segments,  and  a  balance  beam  or  seesaw  that 
provides  a  way  to  represent  accumulation  of  damage. 

At  the  beginning  of  the  program,  a  horizontal  line 
appears  at  the  lower  left  hand  corner  of  the  border  outline 
on  the  screen.   The  appearance  of  this  line  indicates  the 


Figure  6:  'Damage  Model'  Screen  Display  41 

beginning  of  a  lifecycle.  The  line  will  increase  to  a  random 
length  after  which  a  new  horizontal  line  will  appear  a  random 
vertical  distance  above  the  old  one.  The  new  line  will 
increase  to  a  random  length  then  another  jump  will  occur,  and 
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so  on.  The  length  of  a  line  represents  the  interarrival  time 
between  events.  The  occurrence  of  a  jump  represents  an  event 
that  causes  damage,  and  the  vertical  magnitude  of  the  jump 
represents  the  amount  of  damage  incurred.  The  event  arrivals 
form  a  Poisson  process,  and  the  amounts  of  damage  are 
independent  and  identically  distributed  random  variables, 
independent  of  the  event  arrival  times,  (the  damage  amounts 
are  exponential  in  this  model). 


Figure  7:  'Damage  Model'  Screen  Display  #2 

When  a  jump  occurs,  'weight'  is  added  to  the  left  side  of 
the  seesaw.  The  amount  of  weight  added  corresponds  to  the 
amount  of  damage  incurred.  The  dotted  line  near  the  top  of 
the  screen  border  represents  a  threshold;  a  point  which  when 
reached  will  signal  the  end  of  the  lifecycle.  The  final  jump 
to  this  level  will  add  enough  'weight'  to  the  seesaw  to  cause 
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it  to  swing  from  its  initial  position  of  rest  to  the  opposite 
extreme;  a  failure  has  occurred.  The  program  then  resets  the 
screen  images  and  restarts  with  a  new  lifecycle. 

Again,  there  are  few  keyboard  commands.  There  is  no 
fast/slow  clock  speed  option  or  parameter  change  routine. 
Those  options  available  are  common  to  all  programs  and  will 
be  explained  in  the  concluding  paragraphs  that  follow. 

E.   ADDITIONAL  KEYBOARD  COMMANDS 

In  addition  to  the  keyboard  commands  noted  in  the  program 
descriptions,  there  are  two  that  have  not  been  mentioned 
which  are  common  to  all  the  programs.  The  first  is  the 
option  to  temporarily  stop  a  program.  This  may  be  done  by 
pressing  the  space  bar.  The  program  will  do  nothing  further 
and  the  picture  on  the  screen  will  remain  motionless  until 
the  space  bar  is  pressed  again. 

The  second  and  final  command  is  the  option  to  quit  the 
package  of  programs  altogether.  Pressing  the  'Q'  key  will  do 
it;  the  program  will  end  and  the  screen  will  go  blank.  This 
action  is  effective  in  concluding  a  session  with  this 
package.  Note  however,  that  two  further  steps  should  be 
taken  before  the  computer  is  used  for  further  programming;  1) 
press  the  reset  button  to  disable  the  high  resolution 
character  generator  and  clear  memory  of  the  character  sets 
used  in  the  programs,  and  2)  press  the  key  sequence  'FP'  to 
reset  'himem*.   If  these  two  steps  are  not  carried  out  before 
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the  computer  is  put  to  some  other  use  the  user  may  run  out  of 
memory  sooner  than  expected.  Of  course  an  alternate  and 
perhaps  simpler  approach  to  resetting  the  computer  is  to  turn 
it  off  and  then  back  on. 

The  two  tables  that  follow  provide  of  summary  of  keyboard 
commands.  Table  1,  gives  the  definitions  of  each  command. 
Table  2,    supplies  commands  available  in  each  program. 
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Table  1:  Definition  of  Keyboard  Commands 


Key 


Definition 


A 
B 
C 
D 
D 
F 
N 

N 
0 
Q 
S 
T 
U 
U 
Y 

spc  bar 


Run  Poisson  process  for  auto  only 

Run  Poisson  process  for  bus  only 

Change  parameters 

Turn  statistics  display  on 

Lower  water  level  (Dam  Model) 

Run  program  fast 

1st  time;  No  character  animation 

2nd  time;  no  time  plot 

Normal  program  operation  (Dam  Model) 

Turn  statistics  display  off 

End  session 

Run  program  slow  (default) 

Run  decomposed  Poisson  process  (auto  and  bus) 

Run  superimposed  Poisson  process  (auto  and  bus] 

Raise  water  behind  dam  (Dam  Model) 

Resume  graphics  (after  'N*  is  pressed  one  or 
more  times) 

Stop  program  temporarily/resume  program  after 
temporary  stop 
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Table  2:  Reference  Table  for  Keyboard  Commands 

Program  Keyboard  Commands  Available 

Poisson  Process  A,3 ,C ,D, F,N ,0,Q,S ,T,U , Y ,spc  bar 

MMIQ  C,D,F,N,0,S,Q,Y,spc  bar 

Dam  Model  D,F,N ,S ,Q,U , spc  bar 

Damage  Model  QfSpc  bar 
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APPENDIX  A 


PROGRAM  LISTINGS 


The  program  listings  are  written  in  Apple  II  Applesoft 
BASIC.  The  programs  were  listed  using  XLISTER,  a  Beagle 
Brothers,    Inc.    program. 
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A.   POISSON  PROCESS  PROGRAM 


10  GOSUB  10000 

:  REM  INITIALIZE 

999  REM  MAIN  ROUTINE;  CLOCK 

1000  IF  TM  =  30  OR  TM  =  0  THEN  POKE  37  ,  (  PEEK  (37)  +  1) 
*  :  IF  TD  =  1  THEN  PRINT  C$ 


1010  IF  TD  >  1  AND  PEEK  (37)  >  14  THEN  POKE  37,15 

*  :  PRINT  C$ 

*  :  POKE  37,14 

*  :  GOTO  1030 


1020   IF  PEEK  (37)  >  14  THEN  POKE  37,14 
1030   IF  TD  =  3  THEN  1060 


1040  IF  Tl  =  0  THEN  VT  =  1 

*  :  GOSUB  1500 

*  :  GOTO  1120 

*  :  REM  MOVE  AUTO 

1050   IF  TD  =  2  THEN  1070 

1060   IF  T2  =  0  THEN  VT  =  2 

*  •  GOSUB  1500 
GOTO  1120 
REM  MOVE  BUS 


1070   IF  NN  =  0  THEN  1110 


1080   IF  TM  /  60  =  INT  (TM  /  60)  THEN  PRINT  MIN$ (L) ; 
*  :  GOTO  1110 


1090   IF  TM  /  5  =  INT  (TM  /  5)  THEN  PRINT  DASH$ (L) ; 
*  :  GOTO  1110 


1100   PRINT  DOT$ (L) ; 


1110   FOR  I  =  1  TO  R 
:  NEXT 
:  IF  TD  =  1  THEN  POKE  36, (  PEEK  (36)  +  1) 


1120  IF  TD  <  >  3  THEN  Tl  =  Tl  -  1 
1130  IF  TD  <  >  2  THEN  T2  =  T2  -  1 
1140   TM  =  TM  +  1 
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IF  TM  =  60  THEN  TM  =  0 


1150  IF  PEEK  (  -  16384)  >  127  THEN  GOSUB  3000 

*  :  REM  KEY  PRESSED  TO  STOP  CLOCK 

1160  TT  =  TT  +  1 

:  IF  TT  >  32765  THEN  PRINT  Y$;P$ 

*  :  GOTO  10 


1170   IF  CI  >  0  THEN  Ml  =  TT  /  CI 
*  :  IF  TD  =  3  THEN  Ml  =  0 


1180   IF  C2  >  0  THEN  M2  =  TT  /  C2 
*  :  IF  TD  =  2  THEN  M2  =  0 


1190 

IF  SD  =  1  THEN  CH  =  PEEK  (36) 

* 

CV  =  PEEK  (37) 

• 

PRINT  Y$;A$"2" 

* 

VTAB  18 

* 

HTAB  11 

• 

PRINT  TT 

* 

VTAB  19 

* 

HTAB  35 

* 

PRINT  FOR  5,2;#M1 

* 

VTAB  20 

* 

HTAB  35 

* 

PRINT  FOR  5,2;#M2;A$"1" 

* 

GOSUB  2600 

* 

IF  TD  =  1  THEN  POKE  36 , (  PEEK 

(36)  - 

-  1) 

1200   GOTO  1000 

1499  REM  EVENT  SUBROUTINES 

1500  IF  NN  =  0  THEN  2000 


1510   IF  TD  =  2  OR  TD  =  4  THEN  PRINT  DIAM$(1); 
*  :  GOTO  2000 


1520   IF  TD  =  3  THEN  PRINT  DIAM$(8); 
*  :  GOTO  2000 


1530   K  =  4 

:  IF  VT  =  2  THEN  K  =  1 


1540   IF  TM  /  60  =  INT  (TM  /  60)  THEN  PRINT  DIAM$ (K  +  3); 
*  :  GOTO  2000 


1550   IF  TM  /  5  =  INT  (TM  /  5)  THEN  PRINT  DIAM$ (K  +  2); 
*  :  GOTO  2000 
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1560   PRINT  DIAM$(K  +  1); 

1999   REM    VEHICLE  RUN  SUBROUTINE 


2000   CH  =  PEEK  (36) 
CV  =  PEEK  (37) 
PRINT  Y$ 
REM  SAVE  CURSOR  POSITION 


2010   IF  TD  <  >  3  AND  VT 
*  :  Ml  =  TT  /  CI 


2030 
* 

* 

•k 
* 

* 


1    THEN  CI  =  CI  +  1 


2020   IF  TD  <  >  2  AND  VT  =  2  THEN  C2 
*  :  M2  =  TT  /  C2 


=  C2  +  1 


IF  SD  =  1  THEN  PRINT  A$"2" 

VTAB  19 

HTAB  11 

PRINT  CI 

VTAB  19 

HTAB  3  5 

PRINT  FOR  5,2;#M1 

VTAB  20 

HTAB  11 

PRINT  C2 

VTAB  20 

HTAB  35 

PRINT  FOR  5,2;#M2;A$"1" 


2040   IF  N  =  0  THEN  2500 

2050   FOR  I  =  1  TO  33 
VTAB  3 
HTAB  I 

PRINT  CAR$ (VT) ; 
NEXT 
VTAB  3 
HTAB  3  3 
PRINT  E$ 
VTAB  4 
HTAB  33 
PRINT  S$ 

2499  REM  NEXT  EVENT  SUBROUTINE 

2500  IF  VT  =  2  THEN  2540 

2510   T  =  -  El  *  LOG  (  RND  (1) ) 


2520   Tl  =  INT  (T)  +  1 

:  REM  AUTO  ARRIVAL  TIME 
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2530   IF  VT  =  1  THEN  2560 

2540   T  =  -  E2  *  LOG  (  RND  (1)) 


2550  T2  =  INT  (T)  +  1 

:  REM  BUS  ARRIVAL  TIME 

2560  IF  Tl  <  >  T2  THEN  2600 

2570  T  =  RND  (1) 

2580  IF  T  <  P  THEN  T2  =  T2  +  1 

2590  IF  T  >  P  THEN  Tl  =  Tl  +  1 


2600   VTAB  6 
HTAB  5 
PRINT  V$; 
VTAB  16 
HTAB  30 
PRINT  W$; 
POKE  36, CH 
POKE  37, CV 
IF  TD  =  1  THEN  POKE  36, (  PEEK  (36)  +  1) 


2610  RETURN 

2999  REM  KEYPRESS  SUBROUTINE 

3000  CH  =  PEEK  (36) 
:  CV  =  PEEK  (37) 

3010  IF  PEEK  (  -  16384)  =  211  THEN  R  =  1000 

*  :  GOTO  3250 

*  :  REM  RUN  CLOCK  SLOW 


3020   IF  PEEK  (  -  16384)  =  198  THEN  R  =  1 

*  ;  GOTO  3250 

*  :  REM  RUN  CLOCK  FAST 


3030   IF  PEEK  (  -  16384)  =  206  AND  N  =  0  AND  SD 
0 
*  •  PRINT  P$ 
GOTO  3250 
REM  SHOW  STATS  DISPLAY  ONLY 


=  1  THEN  NN  = 


3040   IF  PEEK  (  -  16384)  =  206  THEN  N  =  0 

*  :  GOTO  3250 

*  :  REM  SHOW  ALL  BUT  CHARACTER  ANIMATION 
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3050  IF  PEEK  (  -  16384)  =  217  AND  NN  =  0  THEN  NN  =  1 

*  :  CH    =    0 

*  :  CV   =    5 

*  :  TM   =    0 

*  :  IF  TD  =  1  THEN  CV  =  2 

*  :  GOTO  3250  REM  IN  DISPLAY  ONLY  MODE;  BRING  BACK  TIMING 

PLOT 


3060       IF    PEEK    (    -    16384)    =    217    THEN   N    =    1 

*  :    GOTO    3250 

*  :  REM  NOT  SHOWING  CHARACTER  ANIMATION; BRING  IT  BACK 


3070  IF  PEEK  (  -  16384)  =  196  THEN  SD  =  1 

*  :  GOSUB  3500 

*  :  GOSUB    2600 

*  :  GOTO    3250 

*  :  REM  SET  DISPLAY  MODE 


3080 

IF  PEEK  (  -  16384)  =  207  AND  NN  =  1  THEN  SD  =  0 

*    : 

PRINT  Y$ 

*    : 

VTAB  18 

*    : 

HTAB  1 

*    . 

PRINT  F$ 

*    : 

GOSUB  2600 

*    : 

GOTO    3250 

*    : 

REM  TURN  OFF  DISPLAY 

3090 

IF    PEEK     (    -    16384)    =    212    THEN    TD   =    1 

*    : 

L  =  2 

*    : 

VT  =  0 

•k      . 

POKE    -    16368,0 

*    : 

POP 

*     . 

PRINT  A$"2" 

*     . 

GOTO  4180  REM  DISPLAY  TWO  SETS  OF  TIMING  MARKS 

3100 

IF  PEEK  (  -  16384)  =  194  THEN  TD  =  3 

*    : 

L  =  1 

*    : 

VT  -  1 

*  : 

POKE    -    16368,0 

*  : 

POP 

*     : 

PRINT  A$"2" 

*  : 

GOTO  4180 

*  : 

REM  DISPLAY  BUS  TIMING  MARKS 

3110  IF  PEEK  (  -  16384)  =  193  THEN  TD  =  2 

*  :  L  =  1 

*  :  VT  =  1 

*  :  POKE    -    16368,0 

*  :  POP 

*  :  PRINT  A$"2" 

*  :  GOTO  4180 
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:  REM  DISPLAY  AUTO  TIMING  MARKS 


3120  IF  PEEK  (  -  16384)  =  213  THEN  TD  =  4 

*  ;  L  =  1 

*  :  VT  =  0 

*  :  POKE  -  16368,0 

*  :  POP 

*  :  PRINT  A$"2" 

*  :  GOTO  4180 
REM   DISPLAY  ONE  LINE  FOR  BOTH  VEHICLES 


* 


3130   IF  PEEK  (  -  16384)  =  195  THEN  POP 

*  :  GOTO  4000 

*  :  REM  CHANGE  PARAMETERS 


3140   IF  PEEK  (  -  16384)  =  205  THEN  PRINT  Z$ 
*  :  PRINT  CHR$  (4); "RUN  MENU" 


3150   IF  PEEK  (  -  16384)  =  200  THEN  PRINT  Z$ 
*  :  PRINT  CHR$  (4); "RUN  POISSON  HELP" 


3160   IF  PEEK  (  -  16384)  =  209  THEN  POKE  -  16368,0 

*  :  PRINT  Y$;P$;Z$ 

*  :  END 


3170   CS  =  PEEK  (115)  +  PEEK  (116)  *  256 


3180   IF  PEEK  (  -  16384)  =  178  THEN  PRINT  Z$ 

*  :  PRINT  CHR$  (4);"BL0AD  QUEUE . SET, A" ,CS 

*  :  PRINT  CHR$  (4); "RUN  MMIQ" 


3190   IF  PEEK  (  -  16384)  =  179  THEN  PRINT  Z$ 

*  :  PRINT  CHR$  (4);"BL0AD  DAM. SET , A" ,CS 

*  :  PRINT  CHR$  (4);"RUN  DAM  MODEL" 


3200   IF  PEEK  (  -  16384)  =  180  THEN  PRINT  Z$ 

*  :  PRINT  CHR$  (4);"BL0AD  DAMAGE . SET, A" ,CS 

*  :  PRINT  CHR$  (4); "RUN  DAMAGE  MODEL" 


3210   IF  PEEK  (  -  16384)  =  181  THEN  PRINT  Z$ 

*  :  PRINT  CHR$  (4);"BL0AD  FISH. SET , A" ,CS 

*  :  PRINT  CHR$  (4); "RUN  FISHPOND" 


3220   IF  PEEK  (  -  16384)  <  >  160  THEN  3250 


3230   IF  PEEK  (  -  16384)  =  160  THEN  POKE  -  1636i 
*  :  REM  CLOCK  STOPPED  TO  PAUSE 


3240   IF  PEEK  (  -  16384)  <  >  160  THEN  GOTO  3240 
*  :  REM  AWAITING  SPC  BAR  TO  RESUME 
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3250   POKE  -  16368,0 
POKE  36, CH 
POKE  37, CV 
RETURN 


3499   REM  DISPLAY  SET-UP  SUBROUTINE 


3500   PRINT  Y$;A$"2" 
VTAB  18 
HTAB  5 

PRINT  S$;"TIME=  "  ;  TT 
VTAB  19 
HTAB  1 

PRINT  S$; "COUNT (  )=  ";C1 
VTAB  19 
HTAB  17 
PRINT  "    ="; 
PRINT  FOR  5,2;#G1 

3510   VTAB  19 
HTAB  31 
PRINT  "    ="; 
PRINT  FOR  5,2;#M1 
VTAB  20 
HTAB  1 

PRINT  S$; "COUNT (  )=  ";C2 
VTAB  20 
HTAB  17 
PRINT  "    ="; 
PRINT  FOR  5,2;#G2 
VTAB  20 
HTAB  31 
PRINT  "    ="; 
PRINT  FOR  5,2; #M2 

3520   PRINT  K$;A$"1" 
VTAB  19 
HTAB  7 
PRINT  "5" 
VTAB  19 
HTAB  18 
PRINT  "45" 
VTAB  19 
HTAB  32 
PRINT  "75" 
VTAB  20 
HTAB  7 
PRINT  "6" 
VTAB  20 
HTAB  13 
PRINT  "45" 
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VTAB  20 
HTAB  32 

PRINT  "ve" 

RETURN 


3999   REM  PARAMETER  CHANGE  SUBROUTINE 


4000   PRINT  A$"2";P$ 
:  POKE  -  16368,0 


4010   VTAB  6 
:  HTAB  1 

:  PRINT  S$; "PRESENT  MEANS;  "; C$ ; S$ ;" INTERARRIVAL  ";S$;"TI 
ME  (AUTO) ";C$;S$; "INTERARRIVAL  ";S$;"TIME  (BUS)" 


4020   VTAB  7 
HTAB  26 

PRINT  FOR  5,2; #G1 
VTAB  8 
HTAB  26 
PRINT  FOR  5,2; #G2 

4030   VTAB  10 
HTAB  1 

PRINT  S$;"NEW  MEANS  ?  "; 
GET  R$ 
IF  R$  =  "Y"  THEN  4060 

4040   IF  R$  =  "N"  THEN  4170 

4050   IF  R$  <  >  "N"  THEN  4030 

4060  VTAB  12 

:  HTAB  1 

:  PRINT  F$;S$;" INTERARRIVAL  ";S$;"TIME  (AUTO)  "; 

:  INPUT  "";R$ 


4070   IF  VAL  (R$)  =  0  THEN  4110 


4080   IF  VAL  (R$)  <  =  1  THEN  GOSUB  4500 
*  :  GOTO  4060 


4090   IF  VAL  (R$)  >  99  THEN  GOSUB  4510 
*  :  GOTO  4060 


4100   Gl  =  VAL  (R$) 


4110   VTAB  12 
HTAB  1 
PRINT  F$;S$; "INTERARRIVAL  ";S$;"TIME  (BUS)"; 

INPUT  "";R$ 
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4120   IF  VAL  (R$)  =  0  THEN  4160 


4130   IF  VAL  (R$)  <  =  1  THEN  GOSUB  4500 
*  :  GOTO  4110 


4140   IF  VAL  (R$)  >  99  THEN  GOSUB  4510 
*  :  GOTO  4110 


4150   G2  =  VAL  (R$) 

4160   VTAB  12 
HTAB  1 
PRINT  F$ 
GOTO  4020 

4170  El  =  -  1  /  LOG  (1  -  1  /  Gl) 
:  E2  =  -  1  /  LOG  (1  -  1  /  G2) 
:  P  =  El  /  (El  +  E2) 


4130 

PRINT  P$;Y$ 

4190 

IF  SD  =  1  THEN  VTAB  18 

* 

HTAB  11 

* 

PRINT  "0      " 

* 

VTAB  19 

* 

HTAB  11 

* 

PRINT  "0      " 

* 

VTAB  19 

* 

HTAB  21 

* 

PRINT  FOR  5,2;#G1 

* 

VTAB  19 

* 

HTAB  35 

* 

PRINT  "  0.00" 

* 

VTAB  20 

* 

HTAB  11 

* 

PRINT  "0 

* 

VTAB  20 

* 

HTAB  21 

* 

PRINT  FOR  5,2;#G2 

* 

VTAB  20 

• 

HTAB  3  5 

* 

PRINT  "  0.00" 

4202 

GOSUB  10210 

• 
• 

PRINT  K$;A$"1" 

4210   IF  TD  =  1  THEN  VTAB  3 

*  :  HTAB  1 

*  :  GOTO  1000 
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4220  VTAB  6 
:  HTAB  1 
:  GOTO  1000 

4500   VTAB  12 
HTAB  1 

PRINT  F$;S$;" (VALUE  MUST  BE  >  1)" 
FOR  I  =  1  TO  2000 
NEXT 
RETURN 


4510   VTAB  12 
HTAB  1 

PRINT  F$;S$; " (INPUT  VALUES  <  100)" 
FOR  I  =  1  TO  2000 
NEXT 
RETURN 


9999   REM  INITIALIZE  SUBROUTINE 


10000  A$  =  CHR$  (1) 

B$  =  CHR$  (2) 

C$  =  CHR$  (3) 

D$  =  CHR$  (4) 

E$  =  CHR$  (5) 

F$  =  CHR$  (6) 

K$  =  CHR$  (11) 

L$  =  CHR$  (12) 

0$  =  CHR$  (15) 

P$  =  CHR$  (16) 

Q$  =  CHR$  (17) 

S$  =  CHR$  (19) 

V$  =  CHR$  (22) 

W$  =  CHR$  (23) 

Y$  =  CHR$  (25) 

Z$  =  CHR$  (26) 


10010  DIM  CAR$ (3) ,DOT$ (3) ,DASH$ (3) ,MIN$ (3) ,DIAM$ (9) 

10020  CAR$(1)  =  B$  +  "   ABC"  +  C$  +  "  DEFGHI"  +  Q$  +  D$ 

13030  CAR$(2)  =  B$  +  "  KLMNMPQ"  +  C$  +  "  RSTUVSW"  +  Q$  +  D$ 

10040  DOT$ (1)  =  "X" 

10050  D0T$(2)  =  B$  +  "X"  +  C$  +  "X"  +  Q$  +  D$ 

10060  DASH$ (1)  =  "Y" 

10070  DASH$(2)  =  B$  +  "Y"  +  C$  +  "Y"  +  Q$  +  D$ 
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10080  MIN$(1)  =  "Z" 

10090  MIN$(2)  =  B$  +  "Z"  +  C$  +  "Z"  +  Q$  +  D$ 

10100  DIAM$(1)  =  "1" 

10110  DIAM$(2)  =  B$  +  "X"  +  C$  +  "2"  +  Q$  +  D$ 

10120  DIAM$(3)  =  B$  +  "Y"  +  C$  +  "2"  +  Q$  +  D$ 

10130  DIAM$(4)  =  B$  +  "Z"  +  C$  +  "2"  +  Q$  +  D$ 

10140  DIAM$(5)  =  B$  +  "1"  +  C$  +  "X"  +  Q$  +  D$ 

10150  DIAM$(6)  =  B$  +  "1"  +  C$  +  "Y"  +  Q$  +  D$ 

10160  DIAM$(7)  =  B$  +  "1"  +  C$  +  "Z"  +  Q$  +  D$ 

10170  DIAM$ (8)  =  "2" 


10180  POKE  -  16368,0 

:  PRINT  Z$;0$;A$;0$;D$;P$;A$"1'' 


10190  HCOLOR=  3 

:  HPLOT  1,32  TO  279,32 

10200  El  =  4.48142012 
E2  =  5.48481495 
Gl  =  5 
G2  =  6 
L  =  1 
N  =  1 
NN  =  1 

P  =  El  /  (El  +  E2) 
R  =  1000 
SD  =  0 
SPEED=  255 
TD  =  4 
VT  =  0 

10210  CI  =  0 
C2  =  0 
Ml  =  0 
M2  =  0 
TM  =  0 
TT  =  0 

10220  GOSUB  2500 
VTAB  6 
HTAB  1 
REM   INITIALIZE  FIRST  EVENT 
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10230    RETURN 
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B.   MMIQ  PROGRAM 


10  GOSUB  10000 

:  REM  INITIALIZE 

999  REM  MAIN  ROUTINE; CLOCK 

1000  IF  TM  =  30  OR  TM  =  0  THEN  POKE  37 , (  PEEK  (37)  +  1) 
*  :  PRINT  C$ 


1010   IF  PEEK  (37)  >  15  THEN  POKE  37,15 


1020   IF  Tl  =  T2  AND  Tl  =  0  THEN  EVT  =  1 

*  :  GOSUB  1500 

*  :  GOTO  1100 


1030  IF  SYS  <  >  0  AND  Tl  =  0  THEN  EVT  =  2 

*  :  GOSUB  1500 

*  :  GOTO  1100 

*  :  REM  A  DEPARTURE  EVENT 


1040  IF  T2  =  0  THEN  EVT  =  3 

*  :  GOSUB  1500 

*  :  GOTO  1100 

*  :  REM  AN  ARRIVAL  EVENT 

1050  IF  NN  =  0  THEN  1090 


1060   IF  TM  /  60  =  INT  (TM  /  60)  THEN  PRINT  MIN$ (L) ; 
*  :  GOTO  1090 


1070   IF  TM  /  5  =  INT  (TM  /  5)  THEN  PRINT  DASH$(L); 
*  :  GOTO  1090 


1080  PRINT  DOTS (L) ; 

1090  FOR  I  =  1  TO  R 

:  NEXT 

:  POKE  36, (  PEEK  (36)  +  1) 

1100  T2  =  T2  -  1 

:  IF  SYS  <  >  0  THEN  Tl  =  Tl  -  1 

*  :  REM  IF  NO  ONE  IS  IN  SYSTEM  THEN  CAN'T  HAVE  A  DEPARTURE 
I 


1110  TM  =  TM  +  1 

:  IF  TM  =  60  THEN  TM  =  0 

1120  IF  PEEK  (  -  16384)  >  127  THEN  GCSU3  3000 

*  :  REM  KEY  PRESSED  TO  STOP  CLOCK 
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1130   TT  =  TT  +  1 
QT  =  QT  +  QL 
AQL  =  QT  /  TT 

IF  QT  >  32765  THEN  PRINT  Y$;P$ 
GOTO  10 


1140  IF  SD  =  1  THEN  CH  =  PEEK  (36) 

*  :  CV  =  PEEK  (37) 

*  :  PRINT  Y$;A$"2" 

*  :  VTAB  19 

*  :  HTAB  7 

*  :  PRINT  TT 

*  :  VTAB  2  0 

*  :  HTAB  3  5 

*  :  PRINT  FOR  5 , 2; #AQL; A$ "1 " 

*  :  GOSUB  2760 

*  :  POKE  36, (  PEEK  (36)  -  1) 


1150  GOTO  1000 

1500  IF  NN  =  0  THEN  2000 

1510  IF  EVT  =  1  THEN  PRINT  DIAM$(1); 

*  :  GOTO  2000 


1520   K  =  EVT 

:  IF  EVT  =  3  THEN  K  =  5 


1530   IF  TM  /  60  =  INT  (TM  /  60)  THEN  PRINT  DIAM$ (K  +  2); 
*  :  GOTO  2000 


1540   IF  TM  /  5  =  INT  (TM  /  5)  THEN  PRINT  DIAM$ (K  +  1); 
*  :  GOTO  2000 


1550   PRINT  DIAM$ (K) ; 

2000   CH  =  PEEK  (36) 
CV  =  PEEK  (37) 
PRINT  Y$ 
IF  EVT  =  2  THEN  2400 

2099   REM  ARRIVAL  EVENT  ROUTINE 


2100   L  =  2 

:  IF  SQL  >  25  AND  EVT  =  1  THEN  SYS  =  SYS  +  1 
*  :  GOTO  2400 


2110   IF  SQL  >  25  THEN  SYS  =  SYS  +  1 
*  :  GOTO  2600 
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2120  IF  N  =  0  THEN  2160 

2130  X  =  1 

:  DIST  =  27  -  SQL 

:  IF  SYS  =  0  THEN  DIST  =  29  -  SQL 


2140   GOSUB  2800 

2150   VTAB  3 

HTAB  (X) 
PRINT  BLNK$; 
VTAB  3 
HTAB  (X  +  1) 
PRINT  MAN$  (4)  ; 

2160   SQL  =  SQL  +  2 
:  SYS  =  SYS  +  1 

2170   IF  EVT  =  3  THEN  2600 

2399  REM  DEPARTURE  EVENT  ROUTINE 

2400  IF  N  =  0  THEN  2430 


2410  X  =  29 

:  DIST  =  38 

2420  GOSUB  2800 

:  VTAB  3 

:  HTAB  37 

:  PRINT  CLR$; 

2430  SYS  =  SYS  -  1 

:  IF  SYS  =  0  THEN  L  =  1 

2440  IF  SYS  <  13  THEN  SQL  =  SQL  -  2 

2450  IF  N  =  0  OR  SQL  =  0  THEN  2600 


2460   VTAB  3 
HTAB  26 
PRINT  3LNK$ 
X  =  26 
DIST  =  29 
GOSUB  2800 
VTAB  3 
HTAB  29 
PRINT  BLNK$ 
VTAB  3 
HTAB  30 
PRINT  MAN$ (4 
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2470   IF  SQL  =  2  THEN  2600 

2480   X  =  23 

2490   FOR  I  =  2  TO  SQL  /  2 
:  IF  X  <  1  THEN  2600 

2500   VTAB  3 

HTAB  (X) 
PRINT  MAN$ (1) ; 
FOR  J  =  1  TO  50 
NEXT 
VTAB  3 
HTAB  (X  +  3) 
PRINT  MAN$ (4) ; 
VTAB  3 
HTAB  (X) 
PRINT  BLNK$; 
X  =  X  -  2 
NEXT 

2600   QL  =  SYS  -  1 

:  IF  SYS  =  0  THEN  QL  =  0 

2610   IF  SD  =  1  THEN  VTAB  20 

*  :  HTAB  19 

*  :  PRINT  A$"2";QL;"      " ; A$ " 1 " 


2699  REM  NEXT  EVENT  SUBROUTINE 

2700  IF  EVT  =  3  THEN  2740 
2710  T  =  -  El  *  LOG  (  RND  (1)) 
2720  Tl  =  INT  (T)  +  1 

2730  IF  EVT  =  2  THEN  2760 

2740  T  =  -  E2  *  LOG  (  RND  (1)) 

2750  T2  =  INT  (T)  -i-  1 


27  60   VTAB  7 
HTAB  5 
PRINT  V$; 
VTAB  17 
HTAB  30 
PRINT  W$; 
POKE  36, CH 
POKE  37, CV 


o  ^ 


:  POKE  36, (  PEEK  (36)  +  1) 
:  RETURN 


2800   1=1 

GOSUB  2820 

1  =  2 

GOSUB  2820 

X  =  X  +  1 

1  =  3 

GOSUB  2820 

IF  X  <  DIST  THEN  2800 


2810   RETURN 


2820   VTAB  3 

HTA3  (X) 
PRINT  MANS (I) ; 
FOR  J  =  1  TO  2 
NEXT 
RETURN 

2999  REM  KEYPRESS  SUBROUTINE 

3000  CH  =  PEEK  (36) 
:  CV  =  PEEK  (37) 

3010   IF  PEEK  (  -  16384)  =  211  THEN  R  =  1000 

*  :  GOTO  3210 

*  :  REM  RUN  CLOCK  SLOW 


3020   IF  PEEK  (  -  16384)  =  198  THEN  R  =  1 

*  :  GOTO  3210 

*  :  REM  RUN  CLOCK  FAST 


3030   IF  PEEK  (  -  16384)  =  206  AND  N  =  0  AND  SD  =  1  THEN  NN 
0 

*  :  PRINT  P$ 

*  ;  GOTO  3210 

*  :  REM  SHOW  STATS  DISPLAY  ONLY 


3040   IF  PEEK  (  -  16384)  =  206  THEN  N  =  0 

*  :  GOTO  3210 

*  :  REM  SHOW  ALL  BUT  CHARACTER  ANIMATION 


3050   IF  PEEK  (  -  16384)  =  217  AND  NN  =  0  THEN  NN  =  1 
*  •  TM  =  0 
CV  =  3 
CH  =  0 
GOTO  3210 
REM  IN  DISPLAY  ONLY  MODE;  BRING  BACK  TIMING  PLOT 
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3060   IF  PEEK  (  -  16384)  =  217  THEN  N  =  1 

*  :  GOTO  3600 

*  :  REM  NOT  SHOWING  CHARACTER  ANIMATION;  BRING  IT  BACK 


3070   IF  PEEK  (  -  16384)  =  196  THEN  SD 
*  •  GOSUB  3500 
GOSUB  2760 
GOTO  3210 
REM  SET  DISPLAY  MODE 


3080 

* 

* 
* 

* 
* 


IF  PEEK  (  -  16384)  =  207  AND  NN  =  1  THEN  SD  =  0 

PRINT  Z$ 

VTAB  18 

HTAB  1 

PRINT  F$;0$;T$;A$"1" 

GOSUB  2760 

GOTO  3210 

REM  TURN  OFF  DISPLAY 


3090   IF  PEEK  (  -  16384)  =  195  THEN  POP 

*  :  GOTO  4000 

*  :  REM  CHANGE  PARAMETERS 


3100  IF  PEEK  (  -  16384)  =  205  THEN  PRINT  Z$;0$;B$ 

*  :  VTAB  4 

*  :  HTAB  32 

*  :  PRINT  CLR$;Z$ 

*  :  PRINT  CHR$  (4); "RUN  MENU" 


3110   IF  PEEK  (  -  16384)  =  200  THEN  PRINT  Z$ 
*  :  PRINT  CHR$  (4), -"RUN  QUEUE  HELP" 


3120   IF  PEEK  (  -  16384)  =  209  THEN  POKE  -  16368,0 
*  •  PRINT  Y$;0$;B$ 
VTAB  4 
HTAB  3  2 

PRINT  CLR$;Z$;P$ 
END 


3130   CS  =  PEEK  (115)  +  PEEK  (116)  *  256 


3140   IF  PEEK  (  -  15384)  =  177  THEN  PRINT  Z$;0$;3$ 
*  •  VTAB  4 
HTAB  3  2 
PRINT  CLR$;Z$ 

PRINT  CHR$  (4);"BL0AD  POISSON . SET , A" , CS 
PRINT  CHR$  (4); "RUN  POISSON  PROCESS" 


3150   IF  PEEK  (  -  16384)  =  179  THEN  PRINT  Z$;0$;3$ 

*  :  VTAB  4 

*  :  HTAB  3  2 
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*  :  PRINT  CLR$;Z$ 

*  :  PRINT  CHR$  (4);"BL0AD  DAM.SET, A" ,CS 

*  :  PRINT  CHR$  (4); "RUN  DAM  MODEL" 


3160  IF  PEEK  (  -  16384)  =  180  THEN  PRINT  Z$;0$;B$ 

*  :  VTAB  4 

*  :  HTAB  32 

*  :  PRINT  CLR$;Z$ 

*  :  PRINT  CHR$  (4);"BL0AD  DAMAGE . SET, A" ,CS 

*  :  PRINT  CHR$  (4); "RUN  DAMAGE  MODEL" 


3170  IF  PEEK  (  -  16384)  =  181  THEN  PRINT  Z$;0$;B$ 

*  :  VTAB  4 

*  :  HTAB  32 

*  :  PRINT  CLR$;Z$ 

*  :  PRINT  CHR$  (4);"BL0AD  FISH . SET , A" ,CS 

*  :  PRINT  CHR$  (4); "RUN  FISHPOND" 


3130   IF  PEEK  (  -  16384)  <  >  160  THEN  3210 


3190   IF  PEEK  (  -  16384)  =  160  THEN  POKE  -  16368,0 
*  :  REM  CLOCK  STOPPED  TO  PAUSE 


3200   IF  PEEK  (  -  16384)  <  >  160  THEN  3200 
*  :  REM  AWAITING  SPC  BAR  TO  RESUME 


3210  POKE  -  16368,0 

:  POKE  36, CH 

:  POKE  37, CV 

:  RETURN 


3499   REM  DISPLAY  SET-UP  SUBROUTINE 


3500   PRINT  Y$ 
VTAB  19 
HTAB  1 

PRINT  A$"2";S$; "TIME=  "  ;  TT 
VTAB  20 
HTAB  15 

PRINT  S$;"Q#=  " ; QL 
VTAB  20 
HTAB  27 

PRINT  S$;"AVG  ";S$;"Q#=  ";K$; 
PRINT  FOR  5,2;#AQL;A$"1" 
RETURN 


3599  REM  CLEANUP  SUBROUTINE 

3600  PRINT  Y$; 

:  FOR  I  =  1  TO  30  STEP  2 

:  VTAB  3 
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HTAB  (I) 
PRINT  BLNK$; 
NEXT 


3610   IF  SQL  =  0  THEN  3640 

3620   VTAB  3 
:  HTAB  30 
:  PRINT  MAN$ (4) 

3630  IF  SQL  >  2  THEN  X  =  26 

*  :  FOR  I  =  4  TO  SQL  STEP  2 

*  :  VTAB  3 

*  :  HTAB  (X) 

*  :  PRINT  MAN$ (4) ; 

*  :  X  =  X  -  2 

*  :  NEXT 

3640   POKE  -  16368,0 
VTAB  7 
HTAB  5 
PRINT  V$; 
VTAB  17 
HTAB  30 
PRINT  W$; 
POKE  36, CH 
POKE  37, CV 
RETURN 

3999   REM  PARAMETER  CHANGE  SUBROUTINE 


4000   PRINT  A$"2";P$ 
:  POKE  -  16368,0 


4010   VTAB  8 
:  HTAB  1 

:  PRINT  S$; "PRESENT  "; S$ ; "MEANS; "; C$ ; S$ ; "SERVICE  ";S$;"T 
IME";C$;S$; " INTERARRIVAL  ";S$; "TIME" 


4020   VTAB  9 
HTAB  19 

PRINT  FOR  5,2;  IGl 
VTAB  10 
HTAB  19 
PRINT  FOR  5,2; #G2 

4030   VTAB  12 
HTAB  1 

PRINT  S$;"NEW  ";S$; "MEANS  ?  ";K$; 
GET  R$ 
IF  R$  =  "Y"  THEN  4060 


4040   IF  R$  =  "N"  THEN  4170 

4050   GOTO  4030 

4060   VTAB  14 
HTAB  1 

PRINT  S$;  "SERVICE  ••;S$;"TIME  "; 
INPUT  "";R$ 


4070   IF  VAL  (R$)  =  0  THEN  4110 


4080   IF  VAL  (R$)  <  =  1  THEN  G0SU3  4500 
*  :  GOTO  4060 


4090   IF  VAL  (R$)  >  99  THEN  GOSUB  4510 
*  :  GOTO  4060 


4100   Gl  =  VAL  (R$) 


4110  VTAB  14 

:  HTAB  1 

:  PRINT  F$;S$; "INTERARRIVAL  ";S$;"TIME  "; 

:  INPUT  "";R$ 


4120   IF  VAL  (R$)  =  0  THEN  4160 


4130   IF  VAL  (R$)  <  =  1  THEN  GOSUB  4500 
*  :  GOTO  4110 


4140   IF  VAL  (R$)  >  99  THEN  GOSUB  4510 
*  :  GOTO  4110 


4150 


2  =  VAL  (R$) 


4160   VTAB  12 
HTAB  1 
PRINT  F$ 
GOTO  4020 

4170   El  =  -  1  /  LOG  (1  -  1  /  Gl) 
:  S2  =  -  1  /  LOG  (1  -  1  /  G2) 


4180  PRINT  P$;Y$ 

:  IF  NN  =  1  THEN  FOR  I  =  1  TO  38  STEP  2 

*  :  VTAB  3 

*  :  HTAB  (I) 

*  :  PRINT  BLNK$ 

*  :  NEXT 


4190   IF  SD  =  1  THEN  VTAB  19 
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* 

•  HTAB  7 

■k 

:  PRINT  "0 

* 

:  VTAB  20 

* 

:  HTAB  19 

* 

:  PRINT  "0    " 

* 

:  VTAB  20 

• 

:  HTAB  35 

* 

:  PRINT  "  0,00" 

4200 

IF  NN  =  1  THEN  PRINT  A$"l" 

* 

.  GOSUB  10250 

* 

:  CH  =  1 

* 

CV  =  7 

* 

:  GOTO  1000 

4210 

GOSUB  10250 
:  CH  =  1 
:  CV  =  7 
:  GOTO  1000 

4500 

VTAB  14 
:  HTAB  1 

:  PRINT  F$;S$;" (VALUE  MUST  BE  >  1)" 

:  FOR  I  =  1  TO  2000 

:  NEXT 

:  VTAB  14 

:  HTAB  1 

:  PRINT  F$ 

:  RETURN 

4510   VTAB  14 
HTAB  1 

PRINT  F$;S$;" (INPUT  VALUES  <  100)" 
FOR  I  =  1  TO  2000 
NEXT 
VTAB  14 
HTAB  1 
PRINT  F$ 
RETURN 


9999   REM  INITIALIZE  SUBROUTINE 


10000 


A$ 

= 

CHR$ 

(1) 

B$ 

= 

CHR$ 

(2) 

c$ 

= 

CHR$ 

(3) 

D$ 

= 

CHR$ 

(4) 

E$ 

= 

CHR$ 

(5) 

F$ 

= 

CHR$ 

(6) 

K$ 

= 

CHR$ 

(11) 

L$ 

= 

CHR$ 

(12) 

0$ 

= 

CHR$ 

(15) 
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P$  =  CHR$  (16) 

Q$  =  CHR$  (17) 

S$  =  CHR$  (19) 

T$  =  CHR$  (20) 

V$  =  CHR$  (22) 

W$  =  CHR$  (23) 

Y$  =  CHR$  (25) 

Z$  =  CHR$  (26) 

10010  DIM  MAN$ (5) ,DOT$ (2) ,DASH$ (2) ,MIN$ (2) ,DIAM$ (7) 

10020  MAN$(1)  =  B$  +  "  A"  +  C$  +  "  BC"  +  C$  +  "  DE"  +  D$ 

10030  MANS  (2)  =  B$  +  "  FG"  +  C$  +  "  HI"  +  C$  -f-  "  JK"  +  D$ 

10040  MAN$(3)  =  B$  +  "LM  "  +  C$  +  "NOP"  +  C$  +  "QR  "  +  D$ 

10050  MANS (4)  =  BS  +  "S"  +  CS  +  "T"  +  CS  +  "U"  +  DS 

10060  MANS (5)  =  BS  +  "22"  +  C$  +  "22"  +  DS 

10070  DOTSd)  =  BS  +  "V"  +  CS  +  "  "  +  QS  +  DS 

10080  DOTS  (2)  =  BS  +  "V"  +  CS  +  "V"  +  QS  +  DS 

10090  DASHS(l)  =  BS  +  "W"  +  CS  +  "  "  +  QS  +  DS 

10100  DASHS(2)  =  BS  +  "W"  +  CS  -1-  "W"  +  QS  +  DS 

10110  MINS(l)  =  BS  +  "X"  +  CS  +  "  "  +  QS  +  DS 

10120  MINS(2)  =  BS  +  "X"  +  CS  +  "X"  +  QS  +  DS 

10130  DIAMS(l)  =  BS  +  "Y"  +  CS  +  "1"  +  QS  +  DS 

10140  DIAMS(2)  =  BS  +  "V"  +  CS  +  "1"  +  QS  +  DS 

10150  DIAMS(3)  =  BS  +  "W"  +  CS  +  "1"  +  Q$  +  DS 

10160  DIAMS(4)  =  BS  +  "X"  +  CS  +  "1"  +  QS  +  D$ 

10170  DIAM$(5)  =  BS  +  "Y"  +  CS  +  "V"  +  Q$  +  DS 

10130  DIAMS(6)  =  BS  +  "Y"  +  CS  +  "W"  +  QS  +  DS 

10190  DIAMS(7)  =  BS  +  "Y"  +  CS  +  "X"  +  QS  +  DS 

10200  CLRS  =  BS  +  "     "  +  CS  +  "     ••  -^  CS  +  "     "  +  DS 

10210  BLNKS  =  BS  +  "    "  +  CS  +  "    "  +  CS  +  "    "  +  D$ 
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10220  POKE  -  16368,0 

10230  PRINT  Z$;0$;B$;A$"1'' 

FOR  I  =  1  TO  38  STEP  2 

VTAB  4 

HTAB  I 

PRINT  CLR$; 

NEXT 

VTAB  4 

HTAB  32 

PRINT  MAN$(5) ; 0$ ; A$ ; 0$ ; D$ ; P$ 

VTAB  4 

HTAB  3  2 

PRINT  MAN$ (5) ; 0$ ; T$ 

HC0L0R=  3 

HPLOT  1,40  TO  279,40 


10240  El  =  4.48142012 
S2  =  5.48481495 
Gl  =  5 
G2  =  6 
N  =  1 
NN  =  1 
R  =  1000 
SPEED=  255 

13250  AQL  =  0 
EVT  =  3 
L  =  1 
QL  =  0 
QT  =  3 
SQL  =  3 
SYS  =  3 
TM  =  3 
TT  =  0 

10260  GOSUB  2730 
VTAB  4 
HTAB  1 
REM  INITIALIZE  FIRST  EVENT 


13  270  RETURN 
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C.   DAM  MODEL  PROGRAM 

10  GOSUB  10000 

:  REM  INITIALIZE 

20  GOTO  2000 

999  REM  WATER  PLOT  ALGORITHM 

1000  Z  =  .159  -  (Y  -  97)  *  .00258 

:  X  =  INT  (279  -  82.8  *  ((  -  LOG  (6.28  *  Z))  "  .5)) 

:  IF  X  =  2  *  (  INT  (X  /  2))  THEN  X  =  X  +  1 


1010  RETURN 

1999  REM  OUTFLOW  PLOT  ROUTINE 

2000  IF  PEEK  (  -  16384)  >  127  THEN  GOSUB  3000 
*  :  REM  KEY  PRESSED  TO  STOP 


2010  HCOLOR=  0 

:  IF  UP  =  1  THEN  2300 

2020  FOR  I  =  0  TO  T2  -  1 

2030  GOSUB  1000 

2040  IF  Y  >  156  THEN  FOR  J  =  1  TO  R 

*  :  NEXT 

*  :  GOTO  2140 


2049  REM  OUTFLOW  SIMULATION 

2050  FOR  J  =  1  TO  R 
:  NEXT 

:  IF  Y  <  =  150  THEN  HPLOT  1,156  TO  34,150 

*  :  HCOLOR=  6 

*  :  HPLOT  1,156  TO  34,150 

*  :  HCOLOR=  0 


2059   REM  RECESSION  OF  OVERFLOW 


2060  IF  Y  <  115  THEN  HPLOT  (45  +  (114  -  Y)  "  2),Y  TO  X,Y 

*  :  HCOLOR=  5 

*  :  HPLOT  X  +  2,Y  TO  279, Y 

*  :  HCOLOR=  0 

*  :  Y  =  Y  +  1 

*  :  LVL  =  0 


2070   IF  Y  =  113  THEN  HPLOT  47,114  TO  37,150 
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*  :  HPLOT  46,114  TO  36,150 

*  :  HPLOT  35,151  TO  2,156 

2080  IF  Y  =  111  THEN  HPLOT  45,114  TO  35,150 

2089  REM  NORMAL  OUTFLOW 


2090  IF  Y  <  130  AND  Y  >  =  115  AND  LVL  =  3  THEN  HPLOT  61, Y  TO 
X  Y 

*  :  HCOLOR=  5 

*  :  HPLOT  X  +  2,Y  TO  279, Y 

*  :  HCOLOR=  0 

*  :  LVL  =  0 

*  :  Y  =  Y  +  1 


2100  IF  Y  <  143  AND  Y  >  =  130  AND  LVL  =  2  THEN  HPLOT  61, Y  TO 
X,Y 

*  :  HCOLOR=  5 

*  :  HPLOT  X  +  2,Y  TO  279, Y 

*  :  HCOLOR=  0 

*  :  LVL  =  0 

*  :  Y  =  Y  +  1 


2110   IF  Y  >  150  THEN  HPLOT  1,Y  TO  34, Y 


2120   IF  Y  >  =  143  THEN  HPLOT  61, Y  TO  X,Y 
*  •  HCOLOR=  5 

HPLOT  X  +  2,Y  TO  279, Y 
HCOLOR=  0 
LVL  =  0 
Y  =  Y  +  1 


2130   LVL  =  LVL  +  1 

2140   NEXT 

2150   IF  DWN  =  1  THEN  2000 


2160   IF  PEEK  (  -  16384)  >  127  THEN  GOSUB  3000 
*  :  REM  KEY  PRESSED  TO  STOP 


2199  REM  INFLOW  PLOT  ROUTINE 

2200  IF  UP  =  1  THEN  2300 

2209   REM  MOVE  CLOUD  IN  AND  RAIN 


2210   FOR  I  =  1  TO  15 
VTAB  1 
HTAB  1 
PRINT  CLOUDS (I) 
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:  NEXT 


2220   FOR  I  =  1  TO  15 
VTAB  1 
HTAB  I 

PRINT  CLOUDS (15) 
NEXT 

2230   IF  Tl  <  5  THEN  2260 


2240 

FOR  I  =  1  TO 
VTAB  5 
HTAB  21 
PRINT  LIGHTS 
VTAB  5 
HTAB  21 
PRINT  LBLNKS 
NEXT 
VTAB  5 
HTAB  21 
PRINT  LIGHTS 

3 

FOR  I  =  1  TO 

800 

NEXT 

VTAB  5 

HTAB  21 

PRINT  LBLNKS 

2250 

IF  Tl  >  15  THEN  FOR  I  = 

=  1  TO  800 

*  ; 

NEXT 

* 

FOR  I  =  1  TO 

3 

* 

VTAB  5 

* 

:  HTAB  21 

* 

:  PRINT  LIGHTS 

* 

:  VTAB  5 

* 

:  HTAB  21 

* 

:  PRINT  LBLNKS 

* 

:  NEXT 

* 

:  VTAB  5 

* 

:  HTAB  21 

* 

:  PRINT  LIGHTS 

* 

:  FOR  I  =  1  TO 

800 

* 

:  NEXT 

* 

:  VTAB  5 

* 

:  HTAB  21 

* 

:  PRINT  LBLNKS 

2260 

FOR  I  =  1  TO 
:  NEXT 
:  VTAB  5 
:  HTAB  16 
:  PRINT  RAINS 

800 
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2299  REM  INFLOW  ROUTINE 

2300  HCOLOR=  6 

2310  FOR  I  =  0  TO  Tl  -  1 

2320  GOSUB  1000 

2330  FOR  J  =  1  TO  R 

:  NEXT 

2340  IF  Y  <  108  THEN  2430 

2349  REM  OVERFLOW  SIMULATION 

2350  IF  Y  <  115  THEN  HPLOT  (45  +  (114  -  Y)  "2) ,Y  TO  X,Y 

*  :  HCOLOR=  5 

*  :  HPLOT  X  +  2,Y  TO  279, Y 

*  :  HCOLOR=  6 

*  :  Y  =  Y  -  1 

*  :  LVL  =  0 


2360   IF  Y  =  113  THEN  HPLOT  47,114  TO  37,150 

*  :  HPLOT  46,114  TO  36,150 

*  :  HPLOT  35,151  TO  2,156 


2370   IF  Y  =  111  THEN  HPLOT  45,114  TO  35,150 
2379   REM  NORMAL  INFLOW  SIMULATION 


2380   IF  Y  <  130  AND  Y  >  =  115  AND  LVL  =  3  THEN  HPLOT  51, Y  TO 
X,Y 
HCOLOR=  5 

HPLOT  X  +  2,Y  TO  279, Y 
HCOLOR=  6 
LVL  =  0 
Y  =  Y  -  1 


2390  IF  Y  <  143  AND  Y  >  =  130  AND  LVL  =  2  THEN  HPLOT  61, Y  TO 
X,Y 

*  :  HCOLOR=  5 

*  :  HPLOT  X  +  2,Y  TO  279, Y 

*  :  HCOLOR=  6 

*  :  LVL  =  0 

*  :  Y  =  Y  -  1 


2400   IF  Y  >  150  THEN  HPLOT  (157  -  Y)  *  5,Y  TO  34, Y 


2410   IF  Y  >  =  143  THEN  HPLOT  61, Y  TO  X,Y 
*  :  HCOLOR=  5 
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*  :  HPLOT  X  +  2,Y  TO  279, Y 

*  :  HCOLOR=  6 

*  :  LVL  =  0 

*  :  Y  =  Y  -  1 

2420   LVL  =  LVL  +  1 

2430   NEXT 

2440   IF  UP  =  1  THEN  2000 

2500   FOR  I  =  5  TO  7 
VTA3  I 
HTAB  16 
PRINT  RBLNK$ 
FOR  J  =  1  TO  100 
NEXT 
NEXT 

2510   FOR  I  =  15  TO  26 
VTAB  1 
HTAB  I 

PRINT  CLOUDS (15) 
NEXT 

2520   FOR  I  =  27  TO  40 
VTAB  1 
HTAB  I 

PRINT  CLOUDS (I  -  11) 
NEXT 

2799  REM  SET  UP  FOR  NEXT  INFLOW 

2800  T  =  -  5  *  LOG  (  RND  (1) ) 
2810   Tl  =  INT  (T)  +  1 

2820  T  =  -  5  *  LOG  (  RND  (1) ) 

2830  T2  =  INT  (T)  +  1 

2840  IF  INIT  =  1  THEN  RETURN 

2850  GOTO  2000 

2999  REM  CLOCK  STOPPED,  FIND  OUT  WHY 


3000   IF  PEEK  (  -  16384)  =  211  THEN  R  =  1000 

*  :  GOTO  3160 

*  :  REM  RUN  CLOCK  SLOW 
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3010   IF  PEEK  (  -  16384)  =  198  THEN  R  =  1 

*  :  GOTO  3160 

*  :  REM  RUN  CLOCK  FAST 


3020  IF  PEEK  (  -  16384)  =  213  THEN  UP  =  1 

*  :  DWN  =  0 

*  ;  GOTO  3160 

*  :  REM  WATER  INFLOW  ONLY 


3030  IF  PEEK  (  -  16384)  =  196  THEN  UP  =  0 

*  :  DWN  =  1 

*  :  GOTO  3160 

*  :  REM  WATER  OUTLOW  ONLY 


3040  IF  PEEK  (  -  16384)  =  206  THEN  UP  =  0 

*  :  DWN  =  0 

*  :  GOTO  3160 

*  :  REM  RESUME  NORMAL  IN/OUT  FLOW 


3050   IF  PEEK  (  -  16384)  =  205  THEN  PRINT  Z$ 
*  :  PRINT  CHR$  (4); "RUN  MENU" 


3060   IF  PEEK  (  -  16384)  =  200  THEN  PRINT  Z$ 
*  :  PRINT  CHRS  (4); "RUN  DAM  HELP" 


3070   IF  PEEK  (  -  16384)  =  209  THEN  POKE  -  16368,0 

*  :  PRINT  Y$;P$;Z$ 

*  :  END 


3080   CS  =  PEEK  (115)  +  PEEK  (116)  *  256 


3090   IF  PEEK  (  -  16384)  =  177  THEN  PRINT  Z$ 

*  :  PRINT  CHR$  (4);"BL0AD  POISSON . SET , A" ,CS 

*  :  PRINT  CHR$  (4); "RUN  POISSON  PROCESS" 


3100   IF  PEEK  (  -  16384)  =  178  THEN  PRINT  Z$ 

*  :  PRINT  CHR$  (4);"BL0AD  QUEUE . SET, A" , CS 

*  :  PRINT  CHRS  (4); "RUN  MMIQ" 


3110   IF  PEEK  (  -  16384)  =  180  THEN  PRINT  Z$ 

*  :  PRINT  CHR$  (4);"3L0AD  DAMAGE . SET , A" ,CS 

*  :  PRINT  CHR$  (4); "RUN  DAMAGE  MODEL" 


3120   IF  PEEK  (  -  16384)  =  181  THEN  PRINT  Z$ 

*  :  PRINT  CHR$  (4);"BL0AD  FISH. SET , A" ,CS 

*  :  PRINT  CHR$  (4); "RUN  FISHPOND" 


3130   IF  PEEK  (  -  16384)  <  >  160  THEN  3160 


3140   IF  PEEK  (  -  16384)  =  160  THEN  POKE  -  16368,0 
*  :  REM  CLOCK  STOPPED  TO  PAUSE 


3150   IF  PEEK  (  -  16384)  <  >  160  THEN  GOTO  3150 
*  :  REM  AWAITING  SPC  BAR  TO  RESUME 


3160   POKE  -  16368,0 
:  RETURN 

9999   REM  INITIALIZE  SUBROUTINE 


10000  A$ 

=  CHR$ 

(1) 

:  B$ 

=  CHR$ 

(2) 

:  C$ 

=  CHR$ 

(3) 

:  D$ 

=  CHR$ 

(4) 

:  K$ 

=  CHR$ 

(11) 

:  L$ 

=  CHR$ 

(12) 

1  P$ 

=  CHR$ 

(16) 

QS 

=  CHR$ 

(17) 

s$ 

=  CHR$ 

(19) 

Y$ 

=  CHR$ 

(25) 

z$ 

=  CHR$ 

(26) 

10010  DIM  CLOUD$ (29) 

10020  CLOUDS  (1)  =  B$  +  "  "  +  C$  +  "  "  +  C$  +  L$  +  "A"  +  C$  + 
"E"  +  K$  +  D$ 


10030  CL0UD$(2)  =  3$  +  "  "  +  C$  +  "Q"  +  C$  +  "Z"  +  L$  +  "A"  + 
C$  +  "CE"  +  K$  +  D$ 


10040  CLOUDS (3)  =  B$  +  "  "  +  C$  +  "?Q"  +  C$  +  "  Z"  +  L$  +  "A" 
+  C$  +  "DCE"  +  K$  +  D$ 


10050  CLOUDS (4)  =  3S  +  "  "  +  CS  +  "OPQ"  +  CS  +  "Y  Z"  +  L$  +  " 
A"  +  CS  +  "DDCE"  +  KS  +  D$ 


10060  CLOUDS  (5)  =  BS  +  "H"  +  CS  +  "NOPQ"  +  C$  +  "XY  Z"  +  L$  + 
"A"  +  CS  +  "DDDCE"  +  KS  +  DS 


10070  CLOUDS  (6)  =  3S  +  "GH"  +  CS  +  "  NOPQ"  +  C$  +  "WXY  Z"  +  L 
S  +  "A"  +  CS  +  "DDDDCE"  +  KS  +  DS 


10080  CLOUDS (7)  =  BS  +  "FGH"  +  CS  +  "   NOPQ"  +  CS  +  "  WXY  Z" 
+  LS  +  "A"  +  CS  +  "DDDDDCE"  +  KS  +  DS 


10090  CLOUDS (8)  =  BS  +  "EFGH"  +  CS  +  "    NOPQ"  +  CS  ^  "   WXY 
Z"  +  LS  +  "A"  4-  CS  +  "DDDDDDCE"  +  KS  +  DS 


10100  CLOUDS  (9)  =  BS  +  "DEFGH"  +  CS  +  "M    NOPQ"  +  CS  +  "V   W 
XY  Z"  +  LS  +  "A"  +  CS  +  "CDDDDDDCE"  +  KS  +  D$ 


10110  CLOUDS  (10)  =  BS  +  "CDEFGH"  +  CS  +  "LM    NOPQ"  +  CS  +  " 


o  / 


V   WXY  Z"  +  L$  +  "A"  +  C$  +  "CCDDDDDDCE"  +  K$  +  D$ 


10120  CL0UD$(11)  =  B$  +  "BCDEFGH"  +  C$  +  "KLM  NOPQ"  +  C$  + 
"U  V  WXY  Z"  +  L$  +  "A"  +  C$  +  "CCCDDDDDDCE"  +  K$  + 
D$ 


10130  CL0UD$(12)  =  B$  +  "ABCDEFGH"  +  C$  +  "JKLM    NOPQ"  +  C$ 
+  "TU  V   WXY  Z"  +  L$  +  "A"  +  C$  +  "CCCCDDDDDDCE"  + 
K$  +  D$ 


10140  CL0UD$(13)  =  B$  +  "  ABCDEFGH"  +  C$  +  "IJKLM    NOPQ"  +  C 
$  +  "STU  V   WXY  Z"  +  L$  +  "A"  +  C$  +  "CCCCCDDDDDDCE 
"  +  K$  +  D$ 


10150  CL0UD$(14)  =  B$  +  "   ABCDEFGH"  +  C$  +  "  IJKLM    NOPQ"  + 
C$  +  "RSTU  V   WXY  Z"  +  L$  +  "A"  +  C$  +  "BCCCCCDDDDD 
DCE"  +  K$  +  D$ 


10160  CL0UD$(15)  =  B$  +  "    ABCDEFGH"  +  C$  +  "   IJKLM    NOPQ" 
+  C$  +  "  RSTU  V   WXY  Z"  +  L$  +  "A"  +  C$  +  "  BCCCCCD 
DDDDDCE"  +  K$  +  D$ 


10170  CL0UD$(16)  =  B$  +  "    ABCDEFGH"  +  C$  +  "   IJKLM    NOPQ" 
+  C$  +  "  RSTU  V  WXY  Z"  +  L$  +  C$  +  "  BCCCCCDDDDDDC" 
+  K$  +  D$ 


10180  CLOUDS (17)  =  B$  +  "    ABCDEFGH"  +  C$  +  "   IJKLM    NOP" 
+  C$  +  "  RSTU  V   WXY  "  +  L$  +  C$  +  "  BCCCCCDDDDDD" 
+  K$  +  D$ 


10190  CLOUDS (18)  =  B$  +  "    ABCDEFGH"  +  C$  +  "   IJKLM    NO"  + 
C$  +  "  RSTU  V   WXY"  +  L$  +  C$  +  "  BCCCCCDDDDD"  +  K$ 
+  D$ 


10200  CLOUDS  (19)  =  B$  +  "  ABCDEFGH"  +  C$  +  "  IJKLM  N"  + 
C$  +  "  RSTU  V  WX"  +  L$  +  C$  +  "  BCCCCCDDDD"  +  K$  + 
D$ 


10210  CLOUDS  (20)  =  BS  +  "    ABCDEFG"  +  CS  +  "   IJKLM    "  +  CS 
+  "  RSTU  V   W"  +  LS  +  CS  +  "  BCCCCCDDD"  +  KS  +  DS 


10220  RAINS  =  BS  +  LS  +  "OFGMGFPOFGMGFP"  +  CS  +  "QNFPFFSQNFPF 
FS"  +  CS  +  "QQRSRQRQQRSRQR"  +  KS  +  DS 


10230  CLOUDS  (21)  =  BS  +  "    ABCDEF"  +  CS  +  "   IJKLM   "  +  CS  + 
"  RSTU  V   "  +  LS  +  CS  +  "  BCCCCCDD"  +  K$  +  D$ 


10240  CLOUDS (22)  =  BS  +  "    ABODE"  +  CS  +  "   IJKLM  "  +  CS  +  " 
RSTU  V  "  +  LS  +  CS  +  "  BCCCCCD"  +  KS  +  DS 


10250  CLOUDS  (23)  =  3$  +  "    ABCD"  +  CS  +  "   IJKLM"  -  CS  +  "  R 
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STU  V  +  L$  +  C$  +  "  BCCCCC"  +  K$  +  D$ 

10260  CLOUD$(24)  =  B$  +  "    ABC"  +  C$  +  "   IJKL"  +  C$  +  "  RST 
U  "  +  L$  +  C$  +  "  BCCCC"  +  K$  +  D$ 


10270  CLOUDS  (25)  =  B$  +  "    AB"  +  C$  +  "   UK"  +  C$  +  "  RSTU" 
+  L$  +  C$  +  "  BCCC"  +  K$  +  D$ 


10280  CLOUD$(26)  =  B$  +  "    A"  +  C$  +  "   IJ"  +  C$  +  "  RST"  + 
L$  +  C$  +  "  BCC"  +  K$  +  D$ 


10290  CLOUD$(27)  =  B$  +  "    "  +  C$  +  "   I"  +  C$  +  "  RS"  +  L$ 
+  C$  +  "  BC"  +  K$  +  D$ 


10300  CL0UD$(28)  =  B$  +  "   "  +  C$  +  "   "  +  C$  +  "  R"  +  L$  + 
$  +  "  B"  +  K$  +  D$ 


10310  CLOUD$(29)  =  B$  +  "  "  +  C$  +  "  "  +  C$  +  "  "  +  C$  -r  "  " 
+  D$ 


10320  RAIN$  =  B$  +  L$  +  "OFGMGFPOFGMGFP"  +  C$  +  "QNFPFFSQNFPF 
FS"  +  C$  +  "QQRSRQRQQRSRQR"  +  K$  +  D$ 


10330  LIGHT$  =  B$  +  L$  +  "HI"  +  C$  +  "JK"  +  C$  +  "L"  +  K$  +  D 
$ 


10340  RBLNKS  =  B$  +  "  "  +  D$ 

10  350  LBLNK$  =  B$  +  "   "  +  C$  +  "   "  +  C$  +  "  "  +  D$ 


10360  POKE  -  16368,0 

PRINT  Z$ ; 0$ ; A$ ; 0$ ; D$ ; P$ ; A$  /'  2 " 

VTAB  3 

HTAB  1 

PRINT  S$;"ONE  ";S$; "MOMENT  "; S$ ; "PLEASE ..."; K$ 


10370  PRINT  CHR$  (4);"BL0AD  DAM.PIC , A$2000 " 
10380  PRINT  A$"l" 


10390  LVL  =  0 
R  =  750 
SPEED=  255 
Y  =  139 

10400  INIT  =  1 

:  GOSUB  2800 
:  INIT  =  3 


10410  RETURN 
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D.   DAMAGE  MODEL  PROGRAM 


10  GOSUB  10000 

:  REM  INITIALIZE 

999  REM  MAIN  ROUTINE; CLOCK 

1000  IF  T  =  Tl  THEN  Y  =  Y  -  INC 

*  :  IF  Y  <  =  55  THEN 

*  :  GOSUB  2000 

*  :  X  =  14 

*  :  Y  =  151 

*  :  GOTO  1000 


1010   IF  T  =  Tl  AND  Y  >  55  THEN  GOSUB  2000 
*  :  GOTO  1000 


1020  HPLOT  X,Y 

:  X  =  X  +  1 

:  IF  X  >  265  THEN  X  =  14 

1030  IF  PEEK  (  -  16384)  >  127  THEN  GOSUB  3000 

*  :  REM  KEY  PRESSED  TO  STOP  CLOCK 


1040  T  =  T  +  1 

1050  GOTO  1000 

1999  REM  EVENT  ROUTINE 

2000  IF  WY  =  10  THEN  2030 

2009  REM  ADD  WEIGHT 

2010  FOR  I  =  1  TO  INC 
HPLOT  WX,WY 
WX  =  WX  -  2 

IF  WX  =  8  3  THEN  '^  =  WY  -  2 
WX  =  103 
IF  WY  =  9  THEN  2030 


20  20   NEXT 

2030   IF  Y  >  55  THEN  2500 

2199   REM  DAMAGE  THRESHOLD  MOVEMENT 


2200   VTAB  7 

:  HTAB  X  /  7  +  1 
:  PRINT  "7" 
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HCOLOR=  1 

HPLOT  13,48  TO  266,48 

HCOLOR=  5 

FOR  I  =  13  TO  266  STEP  5 

HPLOT  1,55 

NEXT 

HCOLOR=  3 

FOR  I  =  1  TO  12 

J  =  PEEK  (SS) 

NEXT 


2210   VTAB  2 
HTAB  13 

FOR  I  =  1  TO  3 
PRINT  SEESAW$ (I) ; 
FOR  J  =  1  TO  503 
NEXT 

PRINT  BLNK$; 
NEXT 

PRINT  SEESAWS (4) ; 
FOR  I  =  1  TO  3000 
NEXT 
PRINT  BLNK$;SEESAW$ (0) ; 

2499  REM  SET  UP  FOR  NEXT  EVENT 

2500  T  =  -  15  *  LOG  (  RND  (1) ) 
2510   Tl  =  INT  (T)  +  1 

2520  T  =  -  10  *  LOG  (  RND  (1)) 

2530  INC  =  INT  (T)  +  1 

2540  T  =  0 

2550  IF  Y  >  55  THEN  RETURN 


K$ 


2560   VTAB  7 

HTAB  X  /  7  +  1 

PRINT  L$  +  "Z" 

HCCLOR=  1 

HPLOT  13,48  TO  266,48 

HCOLOR=  5 

FOR  I  =  13  TO  266  STEP  5 

HPLOT  1,55 

NEXT 

HCOLOR=  3 

VTAB  8 

HTAB  3 

PRINT  V$; 
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VTAB  19 
HTAB  36 
PRINT  W$; 
PRINT  P$;Y$ 
WX  =  10  3 
WY  =  29 
RETURN 


2999   REM  KEYPRESS  SUBROUTINE 


3000   IF  PEEK  (  -  16384)  =  209  THEN  POKE  -  16368,0 

*  :  PRINT  Y$;P$;Z$ 

*  :  END 


3010   IF  PEEK  (  -  16384)  =  205  THEN  PRINT  Z$ 
*  :  PRINT  CHR$  (4); "RUN  MENU" 


3020   IF  PEEK  (  -  16384)  =  200  THEN  PRINT  Z$ 
*  :  PRINT  CKR$  (4); "RUN  DAMAGE  HELP" 


3030   CS  =  PEEK  (115)  +  PEEK  (116)  *  256 


3040   IF  PEEK  (  -  16384)  =  177  THEN  PRINT  Z$ 

*  :  PRINT  CHR$  (4);"BL0AD  POISSON . SET, A" ,CS 

*  :  PRINT  CHR$  (4); "RUN  POISSON  PROCESS" 


3050   IF  PEEK  (  -  16384)  =  178  THEN  PRINT  Z$ 

*  :  PRINT  CHR$  (4);"BL0AD  QUEUE . SET , A" ,CS 

*  :  PRINT  CHR$  (4), -"RUN  MMIQ" 


3060   IF  PEEK  (  -  16384)  =  179  THEN  PRINT  Z$ 

*  :  PRINT  CHR$  (4);"BL0AD  DAM. SET , A" ,CS 

*  :  PRINT  CHR$  (4); "RUN  DAM  MODEL" 


3070   IF  PEEK  (  -  16384)  =  181  THEN  PRINT  Z$ 

*  :  PRINT  CHR$  (4);"BL0AD  FISH. SET , A" ,CS 

*  :  PRINT  CHR$  (4); "RUN  FISHPOND" 


3080   IF  PEEK  (  -  16384)  <  >  160  THEN  3110 


3090   IF  PEEK  (  -  16384)  =  150  THEN  POKE  -  16368,0 
*  :  REM  CLOCK  STOPPED  TO  PAUSE 


3100   IF  PEEK  (  -  16384)  <  >  160  THEN  3100 
*  :  REM  AWAITING  SPC  BAR  TO  RESUME 


3110   POKE  -  16368,0 
:  RETURN 

9999   REM  INITIALIZE 
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10000  A$ 

=  CHR$  (1) 

:  B$ 

=  CHR$  (2) 

!  C$ 

=  CHR$  (3) 

:  D$ 

=  CHR$  (4) 

:  K$ 

=  CHR$  (11) 

:  L$ 

=  CHR$  (12) 

:  P$ 

=  CHR$  (16) 

.  Q$ 

=  CHR$  (17) 

:  V$ 

=  CHR$  (22) 

:  W$ 

=  CHR$  (23) 

:  Y$ 

=  CHR$  (25) 

:  Z$ 

=  CHR$  (26) 

1001^ 

)  DIM  SEESAWS (4) 

10020  SEESAWS (0)  =  B$  +  "    "  +  C$  +  "    "  +  C$  +  "890 
ABA"  +  C$  +  "  JKLMNOP      ABA"  +  C$  +  " 
RSTUVWX"  +  Q$  +  D$ 


10030  SEESAWS (1)  =  BS  +  "456"  +  CS  +  "GHG"  +  CS  +  "DEF 
ABA"  +  CS  +  "  JKLMNOP      ABA"  +  CS  +  " 
RSTUVWX"  +  QS  +  DS 


3Q 


3Q 


10040  SEESAWS  (2)  =  BS  +  LS  +  "STS"  4-  C$  +  KS  +  "GHG 

"  +  LS  +  "UVU"  +  CS  +  KS  +  "GHG  ABA"  +  CS 

+  "2YZ"  +  LS  +  "ABCWD      NOP"  +  CS  +  "        YGXIJ 
KLM"  +  KS  +  QS  +  DS 


10050  SEESAWS (3)  =  3S  +  "  "  +  LS  +  "STS"  +  CS  +  " 

UVU  "  +  KS  +  "GHG"  +  CS  +  "ABA  GH 

G"  +  CS  +  LS  +  "NOP      DWHBQR"  +  KS  +  "Yl"  +  CS  + 
+  LS  +  "  MLEJFXGY"  +  K$  +  QS  +  DS 


10060  SEESAWS (4)  =  BS  +  "  456"  +  CS  +  " 

GHG"  +  CS  +  "ABA  DEF"  +  CS  +  "ABA 

ONMLKJ"  +  CS  +  "IWCUTSRQ3"  +  QS  +  DS 


10070  BLNKS  =  BS  +  " 
"  +  CS  +  " 
"  +  CS  + 


"  +  CS  +  " 
"  +  CS  +  " 
"  +  QS  -^  DS 


10080  POKE  -  16368,0 

:  PRINT  Z S ; OS ; AS ; OS ; DS ; P$ ;  AS  ; " 1 " 


10090  HCOLOR=  1 

HPLOT  13,152  TO  267,152  TO  267,48  TO  13,48  TO  13,152 

FOR  I  =  55  TO  150  STEP  10 

HPLOT  11,1  TO  13,1 

NEXT 

FOR  I  =  23  TO  266  STEP  20 

HPLOT  1,152  TO  1,154 
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NEXT 

HCOLOR=  5 

FOR  I  =  13  TO  266  STEP  5 

HPLOT  1,55 

NEXT 

HCOLOR=  3 

VTAB  2 

HTAB  13 

PRINT  SEESAW$(0) 


10100  SS  =  -  16336 
WX  =  10  3 
WY  =  29 

X  =  14 
Y  =  151 

10110  GOSUB  2500 

:  REM  INITIALIZE  FIRST  EVENT 


10120  RETURN 
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E.   FISHPOND  LOGO  PROGRAM 


10  GOSUB  10000 

:  REM  INITIALIZE 

999  REM  MAIN  ROUTINE;  CLOCK 

1000  IF  TT  =  Tl  THEN  GOSUB  2000 
*  :  GOTO  1000 


1010   FOR  I  =  1  TO  R 
:  NEXT 


1020   IF  PEEK  (  -  16384)  >  127  THEN  GOSUB  3000 
*  :  REM  KEY  PRESSED  TO  STOP  CLOCK 


1030   TT  =  TT  +  1 
:  GOTO  1000 

1999  REM  FISH  JUMP  SUBROUTINE 

2000  FOR  I  =  1  TO  HT 
HTAB  X 
VTAB  Y 
PRINT  UP$ 

Y  =  Y  -  1 
NEXT 

2010  HTAB  X 
VTAB  Y 
PRINT  BLNK$ 
HTAB  X 
VTAB  Y 
PRINT  STAY$ 
HTAB  X 
VTAB  Y 
PRINT  BLNK$ 

2020   FOR  I  =  1  TO  HT 
HTAB  X 
VTAB  Y 
PRINT  DOWN$ 

Y  =  Y  +  1 
NEXT 

2030   HTAB  X 
VTAB  Y 
PRINT  BLNK$ 
HTAB  X 
VTAB  Y 
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PRINT  SPLSH$ 
HTAB  X 
VTAB  Y 
PRINT  BLNK$ 


2040   X  =  INT  (30  *  RND  (1)  +  4) 

2049  REM  SET  UP  FOR  NEXT  JUMP 

2050  Y  =  INT  (3  *  RND  (1)  +  20) 


2060   HT  =  INT  (  -  10  *  LOG  (  RND  (1))) 
:  IF  HT  >  5  THEN  HT  =  5 


2070  T  =  -  LOG  (1  -  RND  (1) ) 

2080  Tl  =  INT  (T)  +  1 

2090  TT  =  0 

2100  RETURN 

2999  REM  KEYPRESS  SUBROUTINE 

3000  IF  PEEK  (  -  16384)  =  211  THEN  R  =  1000 
*  :  GOTO  3070 


• 


:  REM  RUN  CLOCK  SLOW 


3010   IF  PEEK  (  -  16384)  =  198  THEN  R  =  1 

*  :  GOTO  3070 

*  :  REM  RUN  CLOCK  FAST 


3020  IF  PEEK  (  -  16384)  =  160  THEN  4000 
3030  IF  PEEK  {  -  16384)  =  193  THEN  4200 
3040   IF  PEEK  (  -  16384)  =  201  THEN  4400 


3050  IF  PEEK  (  -  16384)  =  205  THEN  PRINT  Z$; 

*  :  VTAB  7 

*  :  HTAB  1 

*  :  PRINT  V$; 

*  :  VTAB  19 

*  :  HTAB  4  0 

*  :  PRINT  W$;0$;A$;P$;0$;B$;P$;0$;T$;P$;Z$ 

*  :  PRINT  CHR$  (4); "RUN  MENU" 


3060   IF  PEEK  {  -  16334)  =  209  THEN  POKE  -  16368,0 

*  :  PRINT  Z3;P$ 

*  :  END 
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3070   POKE  -  16368,0 
:  RETURN 

3999   REM  PROGRAM  MENU  SUBROUTINE 


4000   PRINT  Z$ 
VTAB  7 
HTAB  1 
PRINT  V$; 
VTAB  19 
HTAB  40 
PRINT  W$ ; 0$ ; A$ ; P$ ; 0$ ; B$ ; P$ ; 0$ ; T$ ; P$ 


4010   VTAB  8 
:  HTAB  35 

:  PRINT  A$"2";S$;"KEY";C$;C$;  TAB (  3) ; S$ ; "ACKNOWLEDGEMEN 
TS";K$;" A" 


40  2  0  PRINT  TAB(  3) ; S$ ;" INTRODUCTION" ;K$ ; " 
. . I " ; C$ ;  TAB (  3 ) ; S$ ; "MENU" ; K$ ; " . . 
M" 


4030   PRINT  Y$;0$;A$;P$;A$"1" 
:  GOTO  3070 

4099   REM  ACKNOWLEDGEMENTS  SUBROUTINE 


4200  PRINT  Z$;P$ 

:  VTAB  3 

:  HTAB  9 

:  PRINT  "-[  '';S$;  "ACKNOWLEDGEMENTS  ]-" 


4210  VTAB  21 
:  HTAB  11 
:  PRINT  "  [SPC]  TO  CONTINUE" 


4220   VTAB  4 
HTAB  1 
PRINT  V$; 
VTAB  20 

HTAB  4  0 
PRINT  W$; 


4230  VTAB  5 
:  HTAB  3 
:  PRINT  S$;"THE  AUTHOR  WISHES  TO  ACKNOWLEDGE" ; C$ ;  TAB (  3 

) ; "THE  USE  OF  THE  FOLLOWING  SOFTWARE" ; C$ ;  TAB (  3);" 

IN  PRODUCING  THIS  PACKAGE;" 


4240   PRINT  C$;C$;  TAB  (  3) ; K$ ;"' APPLESOFT  TOOL  KIT',  ";L$;"3Y 
";S$; "APPLE  ";C$;  TAB (  3) ; S$ ; "COMPUTER ,  ";S$;"INC. 
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4250   PRINT  C$;  TAB (  3) ; K$ ;"' PRINTER  II',  ";L$;"BY  ";S$;"COMP 
UTER  ";S$; "SYSTEMS  ";C$;  TAB (  3) ; S$ ; "DESIGN" 

4260   POKE  -  16368,0 

4270   IF  PEEK  (  -  16384)  <  >  160  THEN  4270 


4280   PRINT  P$ 

:  VTAB  5 

:  HTAB  3 

:  PRINT  S$;"THE  "; S$ ; "APPLESOFT  TOOL  KIT  PROVIDED" ; C$ ; 
TAB (  3); "THE  MEANS  FOR  BUILDING,  DISPLAYING ,"; C$ ; 
TAB(  3); "AND  ANIMATING  CHARACTERS,  PLUS  A";C$;  TAB ( 
3), -"HOST  OF  OTHER  FEATURES  TO  ENHANCE" 


4290   PRINT  TAB(  3); "PROGRAM  DISPLAY ."; C$ ; C$ ;  TAB (  3);S$;"THE 
";S$; "PRINTER  ";K$;"II  "; L$ ; "PACKAGE  PROVIDED" ; C$ ; 
TAB(  3) ; "FORMATTING  OF  NUMERICAL  OUTPUT  FOUND" ;C$; 
TAB(  3); "IN  TWO  PROGRAMS  OF  THIS  PACKAGE." 


4300   POKE  -  16368,0 

4310   IF  PEEK  (  -  16384)  <  >  160  THEN  4310 


4320   PRINT  Y$;P$ 
:  GOTO  4000 

4399   REM  INTRODUCTION  SUBROUTINE 


4400   PRINT  Z$;P$ 
VTAB  3 
HTAB  11 
PRINT  "-[  ";S$; "INTRODUCTION  ]-" 


4410  VTAB  21 
:  HTAB  11 
:  PRINT  " [SPC]  TO  CONTINUE" 


44  20   VTAB  4 
HTAB  1 
PRINT  V$; 
VTAB  2  0 
HTAB  4  0 
PRINT  W$; 


44  30   VTAB  5 
•  H'T'AB  3 

:  PRINT  SS;"THE  FOUR  PROGRAMS  OF  THIS  PACKAGE" ;C$;  TAB ( 
3); "SIMULATE  SITUATIONS  WHICH  ARE  BASED" ;C$;  TAB (  3 
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);"0N  AN  UNDERLYING  " ; S$ ; "POISSON  PROCESS." 


4440   PRINT  C$;  TAB (  3);S$;"THE  FIRST,  LABELED  '"; S$ ; "POISSON 
" ;S$; -PROCESS' ,";C$;  TAB (  3); "IS  A  SIMPLE  DEMONSTR 
ATION  OF  THE";C$;  TAB (  3 ); "OCCURRENCE  OF  EVENTS  IN 
A  ";S$; "POISSON" ;C$;  TAB (  3); "MANNER.   ";S$;"THE  SE 
COND,  ' ";K$;"MM1Q' ," 


4450  PRINT  TAB(  3)  ; L$ ; "SIMULATES  A  SIMPLE  SINGLE-" ;C$;  TAB  ( 
3); "SERVER  QUEUE  WITH  EXPONENTIAL" ; C$ ;  TAB (  3); "INT 
ERARRIVAL  TIMES  AND  EXPONENTIAL" ; C$ ;  TAB (  3);"SERVI 
CE  TIMES,  {THE";K$;"  M/M/1  "; L$ ; "QUEUE) ." ; 


4460   POKE  -  16368,0 

4470   IF  PEEK  (  -  16384)  <  >  160  THEN  4470 


4480   PRINT  P$ 

:  VTAB  5 

:  HTAB  3 

:  PRINT  S$;"THE  THIRD,  THE  •";S$;"DAM  "; S$ ; "MODEL ' ,  PORT 
RAYS";C$;  TAB  (  3); "ONE  OF  THE  CLASSIC  APPLICATIONS 
OF";C$;  TAB  (  3); "THE  COMPOUND  "; S$ ; "POISSON  PROCESS 
.   ";S$;"THE" 


4490   PRINT  TAB(  3); "FOURTH  PROGRAM,  THE  '"; S$ ; "DAMAGE  ";S$;" 
MODEL'  ,";C$;  TAB  (  3)  /"DEPICTS  THE  COMPOUND  ";S$;"PO 
ISSON  PROCESS";C$;  TAB (  3) /"APPLIED  TO  THE  CUMULATI 
VE  DAMAGE" ;C$;  TAB (  3); "MODEL." 


4500   PRINT 

:  PRINT  TAB(  3);S$;"THS  PROGRAMS  ARE  DESIGNED  TO  PROVIDE 
";C$;  TAB(  3); "A  REAL  TIME  PERCEPTION  OF  ";S$;"POIS 
SON  ";C$;  TAB (  3) ; "PROCESSES  IN  ACTION." 


4510   POKE  -  16368,0 

4520   IF  PEEK  (  -  16384)  <  >  160  THEN  4520 


4530   PRINT  P$ 
VTAB  5 
HTAB  3 
PRINT  S$;"THE  MENU  PROVIDES  THE  DETAILS  OF";C$;  TAB (  3 

);"HOW  TO  RUN  EACH  OF  THE  PROGRAMS ."; C$ ;  TAB (  3);S$ 

;"ANY  PROGRAM  MAY  BE  RUN  FROM  ANY" 


4543   PRINT  TAB (  3); "OTHER  PROGRAM  BY  PRESSING  THE  KEY";C$; 
TAB(  3) ; "DESIGNATED  IN  THE  MENU." 


4550   PRINT  C$;  TAB (  3) ; S$ ;" FINALLY ,  EACH  PROGRAM  IS  PROVIDED 
";C$;  TAB(  3); "WITH  ITS  OWN  '";S$;"HELP'  PROGRAM  WH 
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ICH";C$;  TAB (  3) /"DEFINES  KEYBOARD  COMMANDS.   ";S$; 
" PRESS ";€$;  TAB (  3), -"THE  "  ; 


4560   PRINT  ""';S$;"H'  KEY  TO  GET  '  "  ;  S$  ;  "HELP  .  '  " 

4570   POKE  -  16363,0 

4580   IF  PEEK  (  -  16384)  <  >  160  THEN  4580 


4590   PRINT  Y$;P$ 
:  GOTO  4000 

4699   REM  TITLE  SUBROUTINE 

47  0  0   PRINT  0$;T$;P$;0$;B$;A$"1" 
VTAB  4 
HTAB  3 

PRINT  ID$;C$;C$;C$;  TAB (  11 ) ; L$ ; A$ " 2 " ; "BY 
";S$; " DAVISON ";C$;C$ ; 


";K$;"R.  J. 


4710   PRINT  C$;  TAB (  3) ; S$ ; "SUBMITTED  IN  PARTIAL  FULFILLMENT" 
;C$;  TAB(  3); "OF  THE  REQUIREMENTS  FOR  THE  DEGREE" ;C 
$;  TAB(  3);  "OF  ";S$; "MASTER  OF  "; S$ ; "SCIENCE  IN  ";S 
$; "OPERATIONS" ;C$;  TAB (  3) ; S$ ; "RESEARCH  FROM  THE  "; 


4720   PRINT  S$; "NAVAL  "; S$ ; "POSTGRADUATE" ; C$ ;  TAB (  3);S$;"SCH 
OOL,  ";S$; "MONTEREY,  ";S$;"CA.    "; S$ ; "ADVISORS;  "; 
C$;  TAB(  3);K$;"J.  D.  ";S$;"ESARY  AND  ";K$;"J.  L.  " 
;S$; "ELLIS." 


4730   PRINT  C$;C$;  TAB (  11);"[SPC]  TO  CONTINUE" ; 0$ ;A$ ; P$ ; K$ ; A 
$;"1" 
:  RETURN 


9999   REM  INITIALIZE  SUBROUTINE 


0002 

)  AS 

= 

CHR$  ( 

1) 

B$ 

= 

CHR$  ( 

2) 

c$ 

= 

CHR$ 

3) 

D$ 

= 

CHRS 

4) 

E$ 

= 

CHR$ 

'5) 

K$ 

= 

CHR$ 

.11) 

L$ 

= 

CHR$ 

.12) 

.  0$ 

= 

CHR$ 

:i5) 

•  P$ 

= 

CHR$ 

:i6) 

:  Q$ 

= 

CHRS 

:i7) 

:  S$ 

= 

CHR$ 

:i9) 

:  T$ 

= 

CHR$ 

:20) 

:  V$ 

= 

CHR$ 

[22) 

:  W$ 

= 

CHR$ 

[23) 

:  Y$ 

= 

CHRS 

(25) 
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:  Z$  =  CHR$  (26) 
10010  UP$  =  B$  +  "AB"  +  C$  +  "HI"  +  C$  +  "   "  +  Q$  +  D$ 
10020  STAY$  =  B$  +  "CD"  +  C$  +  "JK"  +  Q$  +  D$ 
10030  DOWN$  =  B$  +  "   "  +  C$  +  "EF"  +  C$  +  "LM"  +  Q$  +  D$ 
10040  SPLSH$  =  B$  +  "    "  +  C$  +  "  NO"  +  C$  +  "PQR"  +  Q$  +  D$ 
10050  BLNK$  =  B$  +  "    "  +  C$  +  "    "  +  C$  +  "    "  +  Q$  +  D$ 


10060    ID$    =    B$    +    "    1    2    4    3    3    2    5         1    6    2    7    %    3    3    %    3 "    -^    C$    + 
"891009     (         3#9$&00&0"    +    Q$4-D$ 


10070  POKE  -  16368,0 

10080  R  =  1000 

:  SPEED=  235 

10090  GOSUB  4700 

:  REM  DISPLAY  TITLE 

10100  GOSUB  2040 

:  REM  INITIALIZE  1ST  JUMP 

10110  RETURN 


F.   MENU  PROGRAM 


1000   A$  =  CHR$ 

(1) 

:  C$  =  CHR$ 

(3) 

:  D$  =  CHR$ 

(4) 

:  K$  =  CHR$ 

(11) 

:  L$  =  CHR$ 

(12) 

:  P$  =  CHR$ 

(16) 

:  S$  =  CHR$ 

(19) 

:  V$  =  CHR$ 

(22) 

:  W$  =  CHR$ 

(23) 

:  Y$  =  CHR$ 

(25) 

:  Z$  =  CHR$ 

(26) 

1010   PRINT  Z$;P$;A$"2" 

.  POKE  -  16368,0 

:  HCOLOR=  3 

:  HPLOT  1,24 

TO  279,24 

HPLOT  1,136  TO  279,136 

VTAB  6 

HTAB  5 

>  PRINT  V$; 

VTAB  16 

HTAB  30 

PRINT  W$; 

2000   VTAB  6 
HTAB  11 

PRINT  "-[  ";K$;"MENU  ]-" 
VTAB  7 
HTAB  26 
PRINT  S$; "PRESS" 


2010  VTAB  3 
:  HTAB  1 
:  PRINT  S$;"POISSON  ";  S$  ;  "PROCESS"  ;  K$  ;  " 1";C$ 

;  "MMIQ 2";C$;S$;  "DAM  ";S$;  "MO 

DEL 3" 


2020   VTAB  11 
:  HTAB  1 

:  PRINT  S$;  "DAMAGE  ";S$;  "MODEL 4";C$;S$;"F 

ISHPOND  ";S$;  "LOGO 5";C$;S$;  "QUIT";K$; 

" Q  " 


20  30  VTAB  15 
:  HTAB  16 
:  PRINT  S$; "SELECTION  -[  ]-"; 


2040   VTAB  15 
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HTAB  28 

GET  R$ 

VTAB  15 

HTAB  28 

PRINT  VAL  (R$) ; 

FOR  I  =  1  TO  1000 

NEXT 


2050   IF  VAL  (R$)  =  0  THEN  POKE  -  15368,0 

*  :  PRINT  Y$;P$;Z$ 

*  :  END 


2060   I  =  VAL  (R$) 

CS  =  PEEK  (115)  +  PEEK  (116)  *  256 

ON  I  GOTO  3000,3010,3020,3030,3040 

2070   GOTO  2040 


3000   PRINT  Z$ 

:  PRINT  CHR$  (4);"BL0AD  POISSON . SET , A" ,CS 
:  PRINT  CHR$  (4); "RUN  POISSON  PROCESS" 


3010   PRINT  Z$ 

:  PRINT  CHR$  (4);"BL0AD  QUEUE .SET , A" ,CS 
:  PRINT  CHR$  (4); "RUN  MMIQ" 


3020   PRINT  Z$ 

:  PRINT  CHR$  (4);"BL0AD  DAM. SET , A" ,CS 
:  PRINT  CHR$  (4); "RUN  DAM  MODEL" 


3030   PRINT  Z$ 

:  PRINT  CHR$  (4);"BL0AD  DAMAGE . SET , A" ,CS 
:  PRINT  CHR$  (4);"RUN  DAMAGE  MODEL" 


3040   PRINT  Z$ 

:  PRINT  CHR$  (4);"BL0AD  FISH. SET , A" ,CS 
:  PRINT  CHR$  (4); "RUN  FISHPOND" 


83 


G.   SAMPLE  HELP  PROGRAM 


1000 


2030 
2040 
2050 


A$ 

= 

CHR$ 

(1) 

c$ 

= 

CHR$ 

(3) 

D$ 

= 

CHR$ 

(4) 

K$ 

= 

CHR$ 

(11) 

L$ 

= 

CHR$ 

(12) 

P$ 

= 

CHR$ 

(16) 

s$ 

= 

CHR$ 

(19) 

V$ 

= 

CHR$ 

(22) 

W$ 

= 

CHR$ 

(23) 

Y$ 

= 

CHR$ 

(25) 

z$ 

= 

CHR$ 

(26) 

1010   PRINT  Z$;P$;A$"2" 
POKE  -  16368,0 
HCOLOR=  3 

HPLOT  1,24  TO  279,24 
HPLOT  1,136  TO  279,136 
VTAB  6 
HTAB  5 
PRINT  V$; 
VTAB  16 
HTAB  3  2 
PRINT  W$; 

2000  VTAB  6 
:  HTAB  9 
:  PRINT  "-[  ";K$;"HELP  ]-" 

2010   VTAB  3 

:  HTAB  1 

:  PRINT  S$;"THIS  PROGRAM  PRESENTS  TWO  INDE-  PENDENT  ";S$ 
;"POISSON  PROCESSES.        "; S$ ;" INITIALLY  THEY  ARE 
COMPOSED  AS   ONE.   ";S$;"THEY  MAY  BE  DECOMPOSED  AND 
RUN  TOGETHER  OR  RUN  SEPARATELY." 


2  0  20   VTAB  15 
:  HTAB  1 
:  PRINT  S$;  " 


[SPC]  TO  CONTINUE"; 
POKE  -  15368,0 
IF  PEEK  (  -  16384)  <  >  160  THEN  2040 


PRINT  P$ 

VTAB  6 

HTAB  3 

PRINT  "-[  ";K$; "KEYBOARD  COMMANDS  ]-" 

VTAB  7 
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:  HTAB  2  5 

:  PRINT  S$;"KEY" 

2060   VTAB  8 

:  HTAB  1 

:  PRINT  S$;  "CLOCK  "  ;  S$  ;  "FAST"  ;  K$  ;  " F";C$;S 

$; "CLOCK  ";S$; "SLOW  (DEFAULT) ";K$; " S";C$;S$;"M 

ENU"  ;  K$  ;  " M"  ;  C$  ;  S$  ;  " QUIT"  ;  K$  ;  "  . 

Q" 


2070   VTAB  12 
:  HTAB  1 

:  PRINT  S$;  "CHANGE  ";  S$  ;  "MEANS"  ;  K$  ;  " C";C$;S 

$;  "PAUSE/"  ;S$;  "RESUME SPC  BAR";C$;C$;  " 

[SPC]  TO  CONTINUE"; 


2080   POKE  -  16368,0 

2090   IF  PEEK  (  -  16384)  <  >  160  THEN  2090 


2100   PRINT  P$ 
VTAB  6 
HTAB  3 

PRINT  "-[  ";K$; "KEYBOARD  COMMANDS  ]-" 
VTAB  7 
HTAB  25 
PRINT  S$; "KEY" 


2110   VTAB  8 

:  HTAB  1 

:  PRINT  S$;"EVENTS  ";  S$ ; "SUPERIMPOSED" ; K$ ; " U";C$;S 

$; "EVENTS  ";S$;  "DECOMPOSED" ; K$ ; " T";C$;S$; "A 

UTO  ";S$;  "EVENTS  ";S$;  "ONLY";K$;  " A";C$;S$; 

"BUS  ";S$;"  EVENTS  "  ;  S$  ;  "ONLY"  ;  K$  ;  " B" 


2120   VTAB  15 
:  HTAB  1 
:  PRINT  S$;  " 


[SPC]  TO  CONTINUE"; 
2130   POKE  -  16368,0 
2140   IF  PEEK  (  -  16384)  <  >  160  THEN  2140 


2150   PRINT  ?$ 
VTAB  6 
HTAB  3 

PRINT  "-[  ";K$; "KEYBOARD  COMMANDS  ]-" 
VTAB  7 
HTAB  26 
PRINT  3$; "KEY" 
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2160  VTAB  8 
:  HTAB  1 
:  PRINT  S$; "STATISTICS  ";  S$  ;  "DISPLAY" ; K$ ; " D";C$; 


S$; "STATISTICS  "; S$ ; "DISPLAY  " 
C$;S$; "VEHICLE  "; S$ ; "DISPLAY  " 
N";C$;S$; "TIME  "; S$ ; "DISPLAY  " 


S$;"OFF";K$;" 0"; 

S$; "OFF";K$; " 

S$; "OFF";K$; " 


2170   VTAB  12 
:  HTAB  1 

:  PRINT  S$; "RESUME  "; S$ ; "VEHICLE  "; S$ ; "DISPLAY" ; K$ ; " 

Y";C$;S$; "RESUME  "S$;"TIME  "; S$ ; "DISPLAY  ";K$;" 

. .Y";C$;C$;L$;"       [SPC]  TO  CONTINUE"; 


2180   POKE  -  16368,0 

2190   IF  PEEK  (  -  16384)  <  >  160  THEN  2190 


2200   PRINT  P$ 
VTAB  6 
HTAB  3 
PRINT  "-[  ";K$; "KEYBOARD  COMMANDS  ]-" 


2210  VTAB  3 
:  HTAB  1 
:  PRINT  S$;"THE  STATISTICS  DISPLAY  MUST";C$;"BE  ON  AND  V 

EHICLE  DISPLAY  OFF" ; C$ ; "BEFORE  THE  TIMING  DISPLAY"; 

C$;"CAN  BE  TURNED  OFF.   "; S$ ; "REVERSE" ; C$ ; "THE  ORDE 

R  TO  RESUME  NORMAL" ; C$ ; "OPERATIONS . " 


2220   POKE  -  16368,0 


2230   VTAB  15 
:  HTAB  1 
:  PRINT  "[SPC]  TO  REVIEW,  ";K$;"[Q]  ";L$;"TO  QUIT"; 


2240   IF  PEEK  (  -  16384)  =  160  THEN  PRINT  P$ 
*  :  POKE  -  16368,0 
GOTO  2000 


* 


2250   IF  PEEK  (  -  16384)  =  209  THEN  POKE  -  16368,0 

*  :  PRINT  Z$ 

*  :  PRINT  CHR$  (4); "RUN  POISSON  PROCESS" 


2260   GOTO  224! 
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H.  STARTUP  PROGRAM 

99  REM  PROGRAM  TO  LOAD  HRCG 

100  HOME 
:  HGR 

110    ONERR  GOTO  370 

120    ADRS  =  0 

130    PRINT  CHR$  (4);"BL0AD  RBOOT" 


140    CALL  520 

:  REM  EXECUTE  RBOOT 


150    ADRS  =  USR  (0),"HRCG" 

199    REM  BRING  IN  HRCG,  ADRES=STARTING  ADDRESS 


200    IF  ADRS  <  =  0  THEN  ADRS  =  ADRS  +  65536 
*  :  REM  MAKE  ADRS  POSITIVE 


210  CS  =  ADRS  -  768  *  2 

:  HIMEM:  CS 

220  CH  =  INT  (CS  /  256) 

:  CL  =  CS  -  256  *  CH 

230  POKE  ADRS  +  7,CL 

:  POKE  ADRS  +  8,CH 


240    PRINT  CHR$  (4);"BL0AD  FISH. SET , A" ; CS 

250    PRINT  CHR$  (4);"BL0AD  ROMAN . SET , A" ; CS  +  768 


300  CALL  ADRS 

:  REM  INITIALIZE  HRCG 

310  PRINT  CHR$  (4);"BL0AD  P2,A$9A07" 

:  REM   LOAD  PRINTER  II 


320    CALL  39431 

:  REM   INITIALIZE  PRINTER  II 


330    POKE  216,0 

:  REM  ONERR  OFF 


340    PRINT  CHR$  (4);"3L0AD  FISHPOND. PIC ,A$4000" 
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350    PRINT  CHR$  (4); "RUN  FISHPOND" 

360    END 

370    PRINT  "UNABLE  TO  LOAD" 


380    POKE  216,0 

:  REM  ONERR  OFF 


390    END 
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